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Abstract. 

In the past few years systems with slow dynamics have attracted considerable 
theoretical and experimental interest. Ageing phenomena are observed during this 
ever-lasting non-equilibrium evolution. A simple instance of such a behaviour is 
provided by the dynamics that takes place when a system is quenched from its high- 
temperature phase to the critical point. The aim of this review is to summarize the 
various numerical and analytical results that have been recently obtained for this case. 
Particular emphasis is put to the field-theoretical methods that can be used to provide 
analytical predictions for the relevant dynamical quantities. 

Fluctuation-dissipation relations are discussed and in particular the concept of 
fluctuation-dissipation ratio (FDR) is reviewed, emphasizing its connection with the 
definition of a possible effective temperature. 

The Renormalization- Group approach to critical dynamics is summarized and the 
scaling forms of the time-dependent non-equilibrium correlation and response functions 
of a generic observable are discussed. From them the universality of the associated FDR 
follows as an amplitude ratio. It is then possible to provide predictions for ageing 
quantities in a variety of different models. In particular the results for Model A, B 
and C dynamics of the O(N) Ginzburg-Landau Hamiltonian, and Model A dynamics 
of the weakly dilute Ising magnet and of the ip 3 theory are reviewed and compared 
with the available numerical results and exact solutions. The effect of a planar surface 
on the ageing behaviour of Model A dynamics is also addressed within the mean-field 
approximation. 
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Nomenclature and notations 

A, B, C,. . . J Dynamical universality classes 

FD Fluctuation-Dissipation 

FDR Fluctuation-Dissipation Ratio 

FDT Fluctuation-Dissipation Theorem 

FT Field Theory 

IR Infra-Red 

MC Monte Carlo (simulation) 

MF Mean Field 

MS Minimal Subtraction 

OT Ordinary Transition (for surface critical behaviour) 

RG Renormalization Group 

RIM Random Ising Model 

T c Critical temperature 

T e g Effective temperature 

SDE Short-Distance Expansion 

SpT Special Transition (for surface critical behaviour) 
t, s Times with t > s, if not differently specified 
TTI Time Translation Invariance 
TRS Time Reversal Symmetry 
UV Ultra- Violet 
(dq) d d q/(2n) d 
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1. Introduction 

Understanding the non-equilibrium dynamics of physical systems is currently one of the 
most challenging problems in statistical physics. Equilibrium statistical mechanics has 
been probably one of the most important achievements during the last century. On 
the other hand, in nature equilibrium is more an exception rather than a rule: real 
systems may persist out-of-equilibrium for several reasons, as external driving forces, 
very slow relaxation etc. Many efforts are currently aiming at achieving a coherent 
theoretical picture of a prototype of such non-equilibrium system: Glassy systems, such 
as structural glasses and spin-glasses. One of the most striking features of glassy systems 
is a dramatic slowing down of relaxational processes. After a perturbation (e.g., a sudden 
temperature change) the slow dynamics manifests itself in one-time quantities, such as 
the energy density, through some slow-decaying dependence on the time elapsed from 
the perturbation. However, even if these quantities get very close to an asymptotic 
value, this does not mean that the system is reaching some stationary or meta-stable 
state. This is elucidated by looking at two-time quantities — mainly correlation and 
response functions — that depend on the two times s and t > s not only through their 
difference, even for long times. Furthermore the decay as a function of t is slower for 
larger s. This phenomenon is usually referred to as ageing [1]: older samples respond 
more slowly. The time s is called the "age" of the system or also waiting time, being 
the time elapsed since the preparation of the system. Useful quantities such as the 
temperature of the system can not be defined in this genuine non-equilibrium regime. 

Ageing behaviour has been experimentally observed in a variety of glassy systems 
(for a review see [2]), and it is expected whenever the equilibration time exceeds any 
experimental observation time-window. This is the case for magnetic systems quenched 
from the high-temperature phase to their critical point or below it, when the critical 
relaxation and the phase ordering take place, respectively. In the latter case both theory 
and experiment clearly display ageing behaviour as reviewed in [3]. 

An important aspect of ageing systems is that the equilibrium fluctuation- dissi- 
pation theorem (FDT) does not hold. Such "violations" of the FDT have been the 
starting point for the introduction of fruitful ideas in the field of glassy systems such 
as non-equilibrium fluctuation-dissipation relations [4] and effective temperatures [5, 6], 
reviewed in [7, 8, 9]. These intriguing theoretical ideas are currently under experimental 
investigation in glassy [10, 11, 12, 13, 14] and granular [15] materials. 

However, as pointed out by Cugliandolo, Kurchan and Parisi [16], "violations" of 
the FDT are not peculiar to glassy and disordered systems: Even a magnetic material 
quenched to its critical point displays non-trivial fluctuation-dissipation (FD) relations. 
Therefore, many aspects of ageing behaviour in classical ferromagnetic models have been 
investigated. Special attention has been devoted to the determination of the so-called 
Fluctuation-Dissipation Ratio (FDR), in particular after the observation [17, 18] that it 
is universal in the sense of Renormalization-Group (RG) theory. Due to universality 
in critical phenomena, this quantity, as well as other aspects of ageing dynamics, 
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is expected to be the same for all the systems belonging to the same universality 
class, irrespective of their specific realization and microscopic details. In particular the 
predictions for universal quantities associated with the dynamics of real systems and 
of the corresponding lattice models can be obtained by means of suitable "mesoscopic" 
field-theoretical models. Understanding ageing behaviour in such simple instances might 
also provide insight into the phenomenology of glassy systems. 

The aim of this paper is to give a comprehensive review of the rather large amount 
of theoretical ideas developed to describe ageing phenomena at the critical point, since 
currently many results are scattered in the literature. Particular emphasis is given to 
the perturbative field-theoretical approach as applied in [19, 20, 21, 22, 23] and to the 
comparison among the results of this analysis, the numerical estimates and the exact 
solutions that have been provided for the simplest models. In this paper we do not 
discuss the vast field of ageing in glassy systems, being well beyond our aims. When 
required, we refer the reader to the comprehensive and updated available reviews [7, 8]. 

The review is organized as follows. 

In section 2, we recall some basic facts concerning the description of dynamics 
in statistical systems. Then, after a simple derivation of the fluctuation-dissipation 
theorem (FDT) in equilibrium, we introduce the FDR in various forms and we remark 
its importance in understanding the dynamics of slow-relaxing systems, and in particular 
its connection with the idea of effective temperature. 

In section 3 we recall the basic steps which allow a path-integral representation 
of the mesoscopic dynamics, given by a suitable Langevin equation. After introducing 
the concept of dynamic universality classes, we briefly summarize the renormalization 
procedure for the corresponding field theories. The scaling forms of the response and 
correlation functions of the order parameter and of other observables are obtained from 
RG equations and short-distance expansion, following the seminal work by Janssen et 
al. [24, 25]. 

In section 4 we consider the purely dissipative dynamics (Model A in the notion 
of [26]) of the Landau- Ginzburg model with O(N) symmetry. We review the numerical 
and exact results available for the relevant observables. We then consider in more detail 
the model within the Gaussian approximation, where the FDRs relative to different 
observables are equal. Finally we review the calculation of two-loop response and 
correlation functions of the order parameter and observables that are quadratic in the 
order parameter. From the associated FDRs it has been shown that the Gaussian 
equivalence of all FDRs does not hold when interactions are taken into account. 

In section 5 we investigate the effects of a spatial surface on the ageing properties of 
Model A dynamics. We report the scaling forms for correlation and response functions 
and we provide a solution of the model in the Gaussian approximation. 

In section 6 the dynamics of a 0(iV)-symmetric conserved order parameter (Model 
B) is analyzed. We review numerical and analytical results. We derive the scaling forms 
for correlation and response functions solving the corresponding RG equations. These 
findings agree with the solutions of particular models and with some recent conjectures. 
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We solve the model exactly in the Gaussian approximation and in the limit of an infinite 
number of components of the order parameter. Beyond the Gaussian approximation we 
show that the FDR is not affected by one-loop corrections. 

In section 7 the case of a non-conserved order parameter coupled to a conserved 
density (Model C) is considered. The results for two-point response and correlation 
functions are provided in a loop expansion up to one loop. 

In section 8 the purely dissipative dynamics (Model A) of a weakly dilute Ising 
model is considered up to the first order in the loop expansion. 

In section 9 the purely dissipative dynamics (Model A) of a y? 3 Landau-Ginzburg 
Hamiltonian is considered up to the first order in e = 6 — d expansion. 

Finally, in section 10, we conclude the review presenting some issues that deserve 
further investigation. 

2. Dynamics, the Fluctuation-Dissipation Theorem and 
Fluctuation-Dissipation Ratios 

2.1. Dynamics 

Let us consider a system in contact with a thermal bath at a given temperature T. In 
principle the dynamics of the system is specified by its microscopic Hamiltonian, either 
classical or quantum, via the evolution equations for the density matrix and phase-space 
density, respectively. However this fully microscopic approach is rarely viable for real 
statistical systems. 

An alternative approach makes use of the representation of the dynamics as a 
stochastic process defined in the space of configurations of the system [27, 28]. At this 
description level, which can still be termed microscopic, one has to assign the transition 
rates between different configurations, on the base of physical considerations. Then the 
master equation rules the time evolution of the probability of finding the system in a 
given state. Denoting by £ the generic configuration of the system and by lf[C m £'] 
the rate of the transition £ i— > (£', the master equation reads 

d t P[€, f ; £ ] = ^{W[(t - C]P[<£, t; £o] - W[€ » <t]P[C, t; <£ ]} , (1) 
e 

where P[€,t; €q] is the probability of finding the system in the configuration € at time 
t starting from the configuration Co at time t — 0. It is possible to show that an 
equilibrium stationary distribution function -P e q[<£] (e.g., the Boltzmann-Gibbs or the 
microcanonical one) is reached for t — > oo if the so-called detailed balance 

W[€ €']P cq [(t] = W[<t £]P cq [£'} , (2) 

holds (a proof of this statement can be found in [27]). Equation (2) has to be fulfilled by 
the transition rates in the master equation in order to describe a dynamics towards an 
equilibrium state. All the dynamical properties of the system can therefore be computed 
once P[£,t; £ ] has been determined by solving the master equation with the assigned 
transition rates. 
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Practically, solving the master equation (1) is rarely feasible. A description of the 
dynamics in terms of mesoscopic variables is in some cases preferable, since it focuses 
directly on those quantities that are expected to determine the dynamical properties 
at length and time scales (referred to as mesoscopic) that are much larger than the 
microscopic ones (at atomic or molecular level) but still small compared to macroscopic 
one set by the dimension of the sample. Mesoscopic variables (or observables) , such 
as the local magnetization density in magnetic systems or the local particle density in 
fluids, are obtained by averaging the corresponding microscopic quantities on mesoscopic 
length and time scales. This average is usually referred to as coarse-graining. 

A viable approach to dynamics, which was first successfully applied to the Brownian 
motion, consists of a description taking advantage of the separation between the typical 
time-scale of fast (microscopic) and slow (mesoscopic) dynamical processes, clearly 
emerging in some cases. It is natural to assume that the dynamics of the mesoscopic 
observables can be described as the result of an effective slow deterministic drift towards 
a stationary state (equilibrium or not) and of a stochastic force that sums up the 
effect of the fast microscopic fluctuations. Of course this description fails to reproduce 
the dynamics taking place at microscopic time and length scales. Hence, when the 
macroscopic physics is crucially related to some microscopic events (as claimed, for 
instance, in the case of rare intermittent events in some glasses [29]) such a mesoscopic 
description is not expected to capture the relevant physical mechanisms. We will strictly 
consider systems for which the mesoscopic description is feasible. 

Let us assume that the mesoscopic properties of the system are described by a set 
of variables fa, labeled by an index %. In terms of {0j} the Hamiltonian is given by Ti.[(f>]. 
The previous heuristic considerations motivate the assumption that the dynamics of the 
system is described by the Langevin equation 

W) = -E D «^+c*(*)- (3) 

The first term in the r.h.s. represents the deterministic evolution and Dy is called kinetic 
coefficient. The stochastic nature of the equation is embodied in the random noise 0(t). 
It is specified by its functional probability distribution function, or equivalently by 
the corresponding moments. In the spirit of the central limit theorem, the first two 
moments are sufficient to characterize the stochastic process, which hence is assumed 
to be Gaussian. Denoting with (•) the mean over the possible realizations of the noise, 
we assume (C?(0) = an d (Ci(t)Cj(^)) = 2Ny<5(t — t'), i.e., a so-called Gaussian white 
noise. 

From the Langevin equation (3) it is possible to derive a Fokker-Planck partial 
differential equation for the probability distribution function P[{fa},t] of fa at time t, 
given an initial condition [30, 31, 28] 

d t P[{fa}A = £r 
5 fa 

The stationary distribution d t P[{fa},t] = is the Boltzmann-Gibbs probability of 



^ T 5P _ 8H „ 
N ifc — + D jfc — P 



(4) 
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equilibrium statistical mechanics Pbg[{0?}] °c e if and only if 

D = /3N. (5) 

This condition is the well-known Einstein's relation. If, given two arbitrary 
configurations {(pi} and {$}, a realization of the noise C« exists such that starting 
from {<f>i} and evolving according to equation (3), the configuration {(p' { } is reached in 
a finite time, the system is said to be ergodic. If ergodicity holds the long-time limit of 
P[{(pi},t] is -Pbg[{0«}] independently of the initial condition. 

The dynamics of the system can be studied by looking at time-dependent quantities 
defined in terms of {(pi} and averaged over the possible realizations of the noise. The 
simplest quantity one might consider among those which are measurable by looking 
at the system only at a given time t (one-time quantities) is ((pi(t)). In the long-time 
limit one expects such a quantity to reach an asymptotic value and therefore it is no 
longer possible to get information about the dynamics of the system. In this sense 
one-time quantities are not usually enough to characterize it. The next step is to 
consider quantities that can be measured by looking at the configuration of the system 
at two different times (two-time quantities). At variance with one-time quantities, 
they also provide information about the dynamics of the system in the long-time limit. 
Interesting two-time quantities are the correlation function of the mesoscopic variables 
{(pi} and the response function. The former is defined as CV,-(t, s) = ((pi(t)(pj(s)) and 
is related to the relaxation of spontaneous (thermal) fluctuations of the variables {(pi} 
whereas the latter expresses the response of the system to an external perturbation. 
In particular let us assume that this perturbation is due to an external field {hi} that 
couples linearly to the variables {(pi} in the Hamiltonian H (i.e., hi is the field conjugate 
to (pi): Hh[(p] = 7~t[<p] — J2i hi(pi. The linear response function Rij(t, s) is defined by 

S((pi(t)) h 



Rij(t,s) 



Shj(s) 



(6) 

h=0 

where (•)/, indicates the average over the stochastic dynamics generated by equation (3) 
with the Hamiltonian Ti h [(j)]. Due to causality, (pi(t) does not depend on hj(s), whenever 
t < s. As a consequence Rij(t, s) oc 9(t — s), where 6{t) — 1 for t > and otherwise. 
In the following we set t > s, so that R(s, t) = 0. 

To calculate the response function, we note that the presence of the external field 
h is equivalent to a shift in the noise 5Q (t) = Dijhj(t), so that equation (6) may be 
written as 

RiAt, s) = V jk (^) = ^DN- 1 )^*)^)) = , (7) 

where we have used Einstein's relation and the algebraic identity (F[(]d(t)) = 
2N;., ( with F[Q an arbitrary function of the noise Ct- A useful relation between 

| This identity is easily derived taking into account that the probability distribution function -P/v[{Ci}] 
of the noise, being Gaussian, satisfies —2NjkSPN/5(k(s) — (j(s)PN- Multiplying both sides of this 
equation by an arbitrary F[Q and performing the functional integration on {Q} (by parts on the l.h.s.) 
yields the result. 
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Rij and CV,- is obtained by considering the time derivative of the correlation function 

(d t - d s )Cij(t,s) = (9 t 0j(O0j(s)) - (0i(t)3,0j(s)) as 

(ft - d s )C t ,(t, s) = -<D ife |^( S )> + (fc^D* J^L) - (0i(OC-(«)> > ( 8 ) 
that can be cast in the form: 

2TR ij (t, s) = (d t - dsldjit, s) - Aij(t, s) , (9) 
where we defined the asymmetry 

4, (i , s) = (*( f )D Jl |M)-(D, l ^( s )). (10) 

Equation (9) relates correlation and response functions by means of the asymmetry. 
A similar equation in terms of microscopic degrees of freedom evolving according to a 
general class of master equations can be found in [32]. 

2.2. The Fluctuation- Dissipation Theorem 

Here we present a simple derivation of the fluctuation-dissipation theorem within the 
approach to dynamics we are discussing. The same result can be obtained in many 
different ways, with different mathematical rigor (see, e.g., [33, 30, 8, 7, 28]). 

When the system equilibrates, equation (9) assumes a simpler form. In fact 
the Time-Translational Invariance (TTI) implies that the correlation and response 
functions satisfy Cij(t,s) = Cij(t — s, 0) and Rij(t,s) = Rij(t — s, 0). Therefore 
(d t — d s )Cij(t,s) = —2d s Cij(t, s). Moreover, equilibrium is characterized by Time- 
Reversal Symmetry (TRS), and so the correlation function of two observables 0\(t) and 
02(0 satisfies (0i(£)0 2 (s)) = (C?i(s)C?2(0)> so that the asymmetry vanishes. Taking 
into account these observations, one concludes that 

R lj (t,s)^f3^^, (11) 

which is just one formulation of the fluctuation-dissipation theorem. Note that the 
FDT holds also if is replaced by the connected correlation function C^\t,s) = 
Cij(t, s) — (4>i(t)){(f)j(s)), since due to TTI (<pi(t)) is independent of time in equilibrium. 

In the following, we are mainly concerned with the case of systems whose physical 
properties close to a critical point are determined by a suitable order parameter (see 
subsection 3.1). In this case 0« and <pj stand for the values of the order parameter 
at different spatial points. Using TTI, space-translational invariance and denoting by x 
the (i-dimensional vector between the two points i and j, equation (11) can be written 
as 

R x (t-s)=/3d s C x (t-s). (12) 

Consider now a local observable 0(x, t) and define its correlation function (always 
with t > s) as C®_ x ,(t, s) = (C?(x, t)0(x', s)) and its response function 

HO(*,t)) h , 



^x-x'(^) S ) 



5ho(x', s) 



(13) 

h o =0 
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where /ie>(x', s) is the field conjugated to 0(x', s). Equation (12) can be generalized to 
O in the form [30, 33] 

R°(t-8)=/3d s C2(t-8), (14) 

i.e., the FDT holds for any observable of the system. This remarkable property allows 
one to define the temperature T = f3~ l through equation (14), independently of the 
specific observable used. 

Whenever the system does not reach thermal equilibrium one can not expect that 
the relations in equations (11), (12), and (14) hold. In this sense their validity may be 
considered as a signal of the fact that the system in equilibrium. Let us consider the 
following experiment. Prepare a system — a glass, a ferromagnet, etc. — in an equilibrium 
state at a high temperature T , greater than the critical or glass transition temperature. 
At time t = quench the system to some temperature T < T by taking it into contact 
with a thermal bath at temperature T and let it evolve. On a general basis, one expects 
that the relaxation towards the equilibrium state corresponding to T is characterized by 
two different regimes: (A) a transient one with non-equilibrium evolution, for t < t eq (T), 
and (B) a stationary regime with equilibrium evolution for t ^> t cq (T), where t eq (T) is 
a characteristic equilibration time of the system. During (A) one expects a dependence 
of the behaviour of the system on initial conditions and both TRS and TTI are broken, 
whereas in (B) TTI and TRS are recovered: The dynamics of fluctuations is given by 
the "equilibrium" one, for which FDT holds. In many systems regime (B) is not reached 
during experimental or even geological times because t eq = oo for all practical purposes. 
These systems always evolve out of equilibrium, even if no perturbation is acting on 
them. In this case the standard concepts of equilibrium statistical mechanics do not 
apply and in particular both TTI and TRS are broken. As a consequence two-time 
quantities, such as the response and correlation functions, depend separately on s and 
t, even for large times. 



2.3. The Fluctuation- Dissipation Ratio 

A measure of the distance from equilibrium of an ageing system evolving in contact with 
a thermal bath at fixed temperature T is the FDR [4] 

^'Sj^jy (15) 

where we assume t > s. When the age s is greater than t eq (T), TRS and TTI hold, and 
the FDT yields A x (t, s) — 1. This is not generically true in the ageing regime. 
The asymptotic value of the FDR 

X°° = lim lim X x=0 (t,s) (16) 

turns out to be a very useful quantity in the description of systems with slow dynamics, 
since X°° = 1 whenever t cq (T) < oo, i.e., when the ageing evolution is interrupted and 
the system crosses over to equilibrium dynamics. Accordingly, X°° ^ 1 is a signal of an 
asymptotic non-equilibrium dynamics. Moreover X°° can be used to define an effective 
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temperature T e g- = T/X°°, which might have some features of the temperature of an 
equilibrium system, e.g., controlling the direction of heat flows and acting as a criterion 
for thermalization [5]. This definition of T c g is closely related to that introduced in 
the context of weak turbulence [34]. Under some assumptions it has been shown that 
X x (t, s) establishes a bridge between the dynamically inaccessible equilibrium state and 
the asymptotic dynamics for large times [35]. 

In equation (12) the FDR has been defined considering the two-time correlation 
function C x (t, s) and the two-time response function R^(t,s) of the order parameter. 
Nevertheless it is possible to define a FDR for any observable O: 

S ) = ' and X ° = lim , lim X *=°^ s). (17) 

O s C^ (t, S) s^oo t^oo 

In [5] it has been shown that a thermometer coupled to the observable O measures, on a 
proper time scale, the temperature T c g = T/Xq, although this point has been source of 
some controversy, since the measured temperature seems to depend on the thermometer 
employed [36]. It has been pointed out (see, e.g., [37]) that the effective temperature 
can be of interest in order to devise some thermodynamics for the system if its value is 
independent of O, i.e., T c g = T cS for all O. This has been explicitly verified for infinite- 
range (mean-field) glass models [5]. Beyond these cases, the observable dependence of 
T e g has been investigated analytically in the case of the trap model [38] and numerically 
for supercooled liquids [39] and in driven systems near a jamming transition [40]. The 
case of ageing in critical systems will be reviewed in section 4. 

As already mentioned, these intriguing features are not only typical of glassy and 
disordered materials, but can be generally found to some extent whenever t eq = oo, i.e., 
in slow-relaxing systems such as ferromagnetic models in the low-temperature phase or 
at the critical point, whose dynamics is characterized by phase ordering and critical 
relaxation, respectively. The phase-ordering dynamics has been mainly investigated 
by means of analytical solutions of exactly solvable models [16, 41, 42, 43, 44, 45] 
and by Monte Carlo simulations of more realistic systems [46, 47, 48, 49, 50]. Due 
to these investigations the phenomenology of phase ordering is quite well understood 
[3], although some open questions remain about the scaling of the response function 
[51, 52, 53, 54, 55, 56]. This issue goes beyond the scope of this review, therefore the 
interested reader is referred to the literature. During phase ordering dynamics, the 
result X°° = (corresponding to an infinite effective temperature, as experimentally 
found in some colloidal glasses [12]) has been put on a firmer ground (see, e.g., [18]). 

Particularly interesting is the case of a quench to the critical point. In fact 
it is possible to get insight into this problem by means of the powerful tools of 
Renormalization Group [30] and scale invariance [57] that have been developed during 
the last decades to investigate mainly equilibrium situations, but that can be applied to 
non-equilibrium relaxation as well. In particular, within this framework, it has been 
argued that X°° and Xq are novel universal quantities of non-equilibrium critical 
dynamics [17, 18, 19, 23]. This clearly emerges from the discussion in sections 3.5.2 
and 3.8. 
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2.4- The Fluctuation- Dissipation Ratio in momentum space 

In the previous section we recalled the definition of the FDR for the response and 
correlation function of a space-dependent observable considered in two points a distance 
x apart, and its long-time limit X°° for x = 0. This is the original form in which the 
FDR has been introduced in the literature on glassy systems [4], mainly with the aim 
of investigating the spatial development of correlations. 

However, within the field-theoretical approach to critical dynamics, it is more 
natural and computationally simpler to focus on the behaviour of observables in 
momentum space. Accordingly hereafter we mainly consider the momentum-dependent 
response i? q (t, s) and correlation C q (t, s) functions, defined as the Fourier transform of 
the corresponding ones in real space either on the lattice or in the continuum. 

To work in momentum space it is worth introducing a quantity that, just like 
X x (t, s), "gauges" the distance from equilibrium evolution. Its natural definition is [19] 

TR^jt, s) 

Note that Xq(t, s) is not the Fourier transform of X x (t, s). The long-time limit 

X£o = lim lim * q=0 (i, s) (19) 

s— >oo t— >oo 

defines, in analogy to X£L , an effective temperature as well. It has been argued that 
these "two" effective temperatures in real and momentum space are equal [19]. Hence, 
the effective temperature may still have a sound physical meaning. Let us outline this 
argument. X x=0 can be written in terms of quantities in momentum space as 

= f(dq)d.C*(t,s) J^^^TR^s) 
x=0 " T/(dg) J R q (t, S ) " J(dq)R n (t,s) — ~ [ « >*, • 1 UJ 

where (dq) = d d q/(2it) d . Thus X~= is the average of Af" 1 , with weight Rq(t, s). In 
the long-time limit t ^> s — > oo, the main contribution to the integral comes from 
the small-q region, provided that Rq(t, s) is peaked around q = with a variance that 
vanishes for t ^> s — > oo. This is the case for most the critical systems we are interested 
in. As a consequence 

XZo = X™ = X 00 . (21) 

To our knowledge this argument is currently the only general one available to support 
equation (21). Several results for specific models [58, 59, 60, 61] confirm equation (21). 
As it stands, the argument is expected to be valid for systems at criticality, whereas the 
extension to other cases might not be straightforward, being fundamental to know the 
form of Rq{t, s) for t ^> s — > oo (see also the discussion at the end of section 8). Given the 
importance of equation (21) for the existence of a proper, single-valued and well-defined 
effective temperature, a more rigorous proof and proper extensions of equation (21) 
would be very welcome. 
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We point out that the definition (18) provides an alternative way to determine 
X°° in numerical simulations, which is expected to increase the accuracy of the results, 
as shown in [58]. The advantage of using X q= o is twofold. On the one hand, 
is defined by using quantities with q = (referred to as coherent observables in [58]) 
which are the sums over the whole sample of the corresponding quantities in real space. 
For instance, when considering the magnetization of a system of volume V, X£L and 
X n=0 are determined by measuring the correlation of the local magnetization m x and 
of the total magnetization M = V^ 1 Xlxev m *> respectively. Even if in the latter case, 
one has to measure connected correlation functions (which are differences and hence 
more noisy than the single terms), one expects that X n= o (thus X^ ) has smaller 
statistical fluctuations than X x=0 (thus X£L ). However the key point is that coherent 
observables directly focus on the non-equilibrated modes of the system, and their FD 
relations typically display very early the crossover to X°°, when the correlation function 
is still comparable to its initial value. In local observables, on the other hand, the various 
relaxing modes are mixed in, and so X°° can only be measured at very large times and 
small ratio s/t, corresponding to quite small values of the correlation function. We will 
show explicitly these differences when discussing the Gaussian model in section 4.2. 

2.5. Quasi- equilibrium regime ofX x (t,s) 

Throughout the review we will be mainly interested in the long-time limit of X x=0 (t, s), 
i.e., in X°°. Nevertheless the function X x=0 (t,s) has also other interesting features in 
different time regimes. 

In particular in the short-time regime, defined as 5t = t — s <C s — > oo, 
general arguments lead to the conclusion that the system is quasi-equilibrated and 
linx^oo X x=0 (s + , s) = 1. This fact can be heuristically accounted for as follows. Let us 
consider the typical quench experiment in which the system is initially in an uncorrelated 
state and then, at t — 0, it is suddenly quenched in the low-temperature phase [3]. 
Thus the domain coarsening takes place. One expects that, at a given time, local 
quantities in two points that are inside the same domain display correlations typical of 
the thermodynamic phase of that domain. Deviations are expected when the two points 
belong to distinct domains. In the latter case the non-equilibrium behaviour is displayed. 
Being the short-time evolution (St <C s) related to processes taking place at short 
spatial distances (at least for systems with short-range interactions), it is equilibrium- 
like, resulting in X x=Q (s + St,s) = 1. Conversely, after a time St comparable with s, 
the behaviour at larger spatial distances starts influencing the dynamics and the system 
falls out of equilibrium. 

When the system is quenched at the critical point no domains are forming. 
Nevertheless, at a given time, the correlation and response functions are equilibrium- 
like when referred to points within a given (growing) range, called dynamic correlation 
length, whereas they depart from equilibrium at larger distances. Thus, even in this 
case, X x=0 (t, s) exhibits a crossover from 1 at t ~ s to X°° for St 3> s. 
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The previous heuristic description can be made more rigorous considering specific 
models. So far quasi-equilibrium behaviour has been observed in essentially all the 
studied cases (see the reviews [7, 8]). 

A general theorem [62], obtainable under quite reasonable assumptions for system 
governed by a dissipative Langevin dynamics, provides a bound for the quantity 
V(t, s) = [1 — X x=0 (t, s)]d s C x=0 (t, s) . This bound implies that V(t, s) — > for fixed and 
finite t — s and s — > oo. Accordingly quasi-equilibrium is reached whenever d s C x=0 (t, s) 
does not vanish. However, for the one-dimensional Ising model, an observable O has 
been found such that X® =0 (s, s) = 3/4 for s — > oo [63], while d s C® =0 (t, s) vanishes in 
the same limit. This puts forward the idea that quasi-equilibrium can be defined only 
on a certain class of observables, termed neutral in reference [63]. 

Nevertheless, the notion of quasi-equilibrium is still very important: It was used 
in [64] to define a nominal thermodynamic temperature for models where the dynamics 
does not satisfy detailed balance and hence T bath , the temperature of the thermal bath, 
does not exist a priori. In any case, this can not be done in general non-equilibrium 
systems that do not satisfy detailed balance. In fact there are several examples of 
models where, due to a different scaling in time of the response and of the correlation 
functions, a non-trivial short-time X x=0 (t, s) can not be defined. This is the case of the 
contact processes [65] and of the deterministic non-linear evolution (i.e., equation (33) 
with C = 0) [66]. 

Let us remark that the quasi-equilibrium regime is not expected to be detectable 
via Xq=o(t,s). Indeed the two-point correlation and response functions involved in its 
definition get non-equilibrium contributions (through the Fourier transform) from the 
corresponding functions in real space taken at points whose distance is bigger than the 
domain size or than the dynamic correlation length. 

3. Field-Theoretical Approach to Non-equilibrium Dynamics 

In this section we briefly review the Renormalization-Group approach to the dynamics 
following a quench from a high-temperature state to the critical point. In particular we 
take advantage of various well-known methods to determine the scaling forms of the non- 
equilibrium response and correlation functions at criticality, computing the associated 
universal quantities in perturbation theory. 

3.1. Critical Phenomena: A brief reminder 

The macroscopic properties of statistical systems (e.g., fluids, magnetic materials, etc.) 
are generally expected to depend strongly on the usually quite complex interactions 
among their microscopic constituents. 

Theoretical insight can be gained only by studying analytically or numerically 
rather simplified models of real systems and the predicted behaviour depends on 
the parameters of the specific model considered. Accordingly, it is important to 
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have a clear understanding of the connection between the model parameters and the 
microscopic properties of the actual system in order to provide predictions comparable 
with experiments. Due to this difficulty only in few cases a quantitative comparison can 
be done. Nevertheless there are circumstances in which a collective behaviour emerges 
which is largely independent of the microscopic details of the actual system and, as a 
consequence, also of the particular model used to describe it. This property is known 
as universality and it naturally characterizes the physical behaviour upon approaching 
a critical point, where the system undergoes a continuous phase transition. 

The onset of a collective behaviour is clearly revealed by the correlation length 
£, defined as the typical distance over which the microscopic variables are correlated. 
Far away from a critical point £ is typically of the order of the range of microscopic 
interactions, whereas it diverges upon approaching the critical point. Accordingly, close 
enough to the transition point, £ becomes mesoscopic and indeed it provides the only 
relevant length-scale of a critical system. 

In view of the universality of critical properties, which is justified within the 
renormalizat ion-group theory (see subsection 3.5.1) and supported by experimental data 
[67], it is possible to study this collective behaviour in terms of suitable mesoscopic 
field-theoretical models, in a formal development of the Landau approach to phase 
transitions [68]. Indeed, as long as one is interested in the behaviour at mesoscopic 
length and time scales, an effective Hamiltonian which reflects the internal symmetries 
of the underlying microscopic system can be used. Such Hamiltonian depends only on 
the order parameter and potentially on a few other slow modes, whose actual nature is 
determined specifically by the system. For instance the order parameter can be identified 
with the magnetization in magnetic materials close to the Curie temperature, with the 
particle density in fluids etc. The previous considerations motivate the assumption that 
the static critical properties of a system with a iV-component vector order parameter 
<f, short-range interactions and O(N) symmetry are captured for large distances by the 
Landau-Ginzburg- Wilson Hamiltonian 



HM= / d d x 



(22) 



where r oc T is a parameter that has to be tuned to the value r 0)C rit in order to approach 
the critical point for T = T c and go > is the coupling constant of the theory. From 
now on we absorb the factor /3 — 1/T in the Hamiltonian, so that 7i henceforth denotes 
the reduced Hamiltonian. Statistical averages of quantities O depending on the order 
parameter field ip are then properly computed, in the equilibrium canonical ensemble, 
as 

«*» ■ 5 -^&s^ ■ 

In a first approximation {mean-field - MF) the integral over the fields is determined 
by the configuration ip MF which minimizes H, i.e., satisfying STi,[(p]/S(p\ v=ipMF = 0. 
Corrections coming from fluctuations around y?MF can be accounted for by considering 



Ageing Properties of Critical Systems 



16 



the successive terms of the formal expansion of 7i[tp] around y?MF- The term which is 
quadratic in ip — <-Pmf defines the so-called Gaussian approximation of Ti. Note that 
the free energy 5 = — In f[d(p]e~ n ^ is given, within the mean-field approximation, 
by Jmf = W[<£>mf]- m this sense TC[(fi\ is sometimes referred to as Landau free-energy 
functional. 

By means of field-theoretical techniques it is possible to determine the non-analytic 
behaviour observed in various thermodynamic quantities and structure factors upon 
approaching the critical point. Such non-analyticities, parameterized by the standard 
critical exponents, some associated amplitude ratios and scaling functions turn out to 
be universal quantities (see, e.g., [67, 69]). The values of universal quantities and scaling 
functions characterize the so-called universality class of the model. 

Upon approaching a critical point the typical time scale of dynamics of the 
fluctuations around the equilibrium state diverges as ~ £ z (critical slowing down), where 
z is the dynamic critical exponent. This provides the natural separation between the 
relevant slow evolution due to the developing collective behaviour and the fast one 
related to microscopic processes. This separation makes the mesoscopic description of 
the dynamics a particularly viable approach to the problem, as explained in section 2.1. 
Indeed it allows one to compute systematically the non-analytic behaviours observed in 
dynamical quantities, e.g., in the low- frequency limit of the dynamic structure factor. In 
turn the associated universal quantities define the dynamic universality class. One finds 
that each static universality class consists of several dynamic sub-universality classes 
which differ, e.g., by different conserved quantities, but nonetheless exhibit the same 
static universal properties. In subsection 3.4 we provide some examples of physical 
systems belonging to the same static universality class but to different dynamical ones. 

3.2. N on- equilibrium critical dynamics 

In the following sections we will focus on a particular instance of non-equilibrium 
behaviour: The one due to a sudden thermal quench to the critical point. From 
the theoretical point of view this behaviour is induced by the initial conditions of the 
dynamics and is not generic. Indeed as soon as the quench is done at a temperature 
slightly above the critical one the system thermalizes in a finite time t eq ~ £, z and reaches 
an equilibrium state characterized by the canonical distribution function proper to the 
mesoscopic Hamiltonian TC of the system. At the critical point the effects of the initial 
conditions persist forever and give rise to a non-equilibrium critical behaviour with some 
universal features. 

It is possible to study theoretically this cases by simply accounting for the initial 
conditions in the mesoscopic evolution equation usually used to describe the dynamics 
of fluctuations around the equilibrium state. In particular we consider only the case 
of a disordered (high-temperature) initial condition. The field-theoretical approach to 
this problem was developed in a seminal paper by Janssen, Schaub and Schmittmann 
[24], who focus on the early stage of the relaxation process after the quench, i.e., on the 
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short-time scaling. The scaling forms obtained there for the two-time correlation and 
response functions C q (t, s) and i? q (i, s) (t > s), in the limit of long times, turn out to 
depend on the ratio s/t. As a consequence, the short-time regime investigated in [24] 
defined by the limit s — > with fixed t is equivalent to t — > oo with fixed s, that is the 
ageing limit of equations (16), (17) and (19). Accordingly we can take advantage of the 
analysis presented in [24] to study ageing behaviour at the critical point. 

The main results obtained from the RG analysis of this problem are the scaling 
forms of the response and correlation functions, in equation (91) and (92), respectively. 
From them one deduces that: 

(i) The critical exponent z appearing in the equilibrium and non-equilibrium scaling 
forms is the same [24], although its physical meaning is different. 

(ii) Only one new independent critical exponent 9 has to be introduced to describe the 
non-equilibrium dynamics [24] , compared to those required for the equilibrium one. 

(iii) The FDR X°° [see equations (15) and (16)] is a universal amplitude-ratio (this 
conclusion was first drawn in [17] on the sole basis of scaling arguments). 

In the next subsections we assume the reader to be familiar with the field-theoretical 
approach to critical phenomena as explained in many introductory textbooks [30, 28, 57] 
and we summarize here only the most relevant points, mainly concerning dynamics. 

3.3. Path- Integral Representation of Dynamics 

The field-theoretical approach to dynamic critical phenomena relies on a path-integral 
description of stochastic processes. Here we recall the basic steps which allow 
the construction of the path-integral associated with a given Langevin equation, as 
developed in [70, 71] and reviewed in [72, 25, 28, 30]. 

Consider the following equation for the field tp (here tp can be either a single field, 
e.g., the order parameter as in Models A and B or a set of slow and conserved modes 
in a more complex case, e.g., tp h- > (tp,e) in Model C - cf. subsection 3.4) 



In this equation T is a local functional of <p(x.,t) and C(x, £) is a zero-mean Gaussian 
white noise with correlation 



where M can either be a constant (as it is the case in Model A) or a differential operator 
acting on x (as in the case of Model B). In all the models we are interested in T has 
the form 



where Ti is a functional playing the role of a reduced Hamiltonian and T> a constant or 
a differential operator. To have TRS and thus a convergence towards the equilibrium 



%>(x,t)=^b(x,*)]+C(x,f). 



(24) 



(C(x,t)C(x',t')) = 2AT5(t-t')5(x-x'), 



(25) 




(26) 
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distribution function for the field P[ip] oc e~ n ^\ the condition V = M has to be fulfilled, 
which is the analogue of equation (5). 

Given an initial condition (pipe, to) the expectation value of a generic observable 
0[ip\ over all possible realizations of the noise ( can be written as 

(O) = J [dC] 0[p c ]P G [Q = J [M 0[ip\ I J [d(]% - ip c )P G [(] } . (27) 

Pg[(] is the Gaussian functional probability distribution function of the noise and ip^ 
is the solution of equation (24) for a given realization of the noise, with the specified 
y?(x, t ) as initial condition. Taking into account that 

8J~~ 

8(ip - <pt) = 5(d t ip - F[ip\ - C) det 



* Sip 



(28) 



it is possible to express the functional 5-function as an exponential by introducing a 
complex auxiliary field ip: 5(ip) = f[dip] exp{ J dtd d x v?(x, t)i[)(x., t)}. Then the average 
over the Gaussian noise is straightforward and leads to 

(O) = Jldipdp]Oe~ St o^ , (29) 

where 

poo p 

S t0 [ip, <p}= dt d d x {<p[d t ip - F[ip\] - ipAfip} . (30) 



Jt J 

The functional S to [ip,ip], usually referred to as the dynamic functional, is the starting 
point for the field-theoretical approach to dynamics. Note that <p(x.,t) has a clear 
physical meaning. Indeed, given an external field h coupled linearly to ip, one has 
Ti[<p, h] = H[ip\ — (3 f d d xhip. This implies, following equation (30), that <p is conjugated 
in S to [ip,<p\ to the external field h, i.e., S to [ip,<p,h] = S to [ip,<p\ — f3 dt f d d x(pVh. As 
a consequence, the linear response of an observable O to the field h is given by 

JM- rm ,, s) VO). (31) 

For this reason <^(x, s) is termed response field. In particular the response function (of 
the order parameter) reads 

S(ip(x,t)) h 



6h(xf, s) 



= P{<p(x!,s)V ( p(x,t)). (32) 



h=0 



From now on we absorb the factor T = in the definition of the response function. 

We remark that in equation (30) the term corresponding to the determinant 
in equation (28) is missing. To be properly evaluated, it requires a discretization 
of the Langevin equation and eventually its expression depends on the chosen 
discretization [28, 73, 25]. Nevertheless the result of the computation of averages such 
as equation (29) is actually independent of the particular choice [72] which can be 
made in such a way to render the determinant equal to one. In turn, this implies that 
R*-x.'{t,t) oc (v?(x, £)C( X ', t)) = 0, corresponding to the so-called Ito prescription in 
stochastic calculus (see, e.g., [27, 74]). Accordingly, in the perturbative expansion of 
averages as (29) all the diagrams with at least one loop of the response function do not 
contribute. 
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3-4- Dynamic Universality Classes 

A classification of several dynamic universality classes was done in the early seventies 
and is reviewed in the classical paper by Hohenberg and Halperin [26]. These universality 
classes have been named with capital letters, from A to J (some new classes have been 
added to the original classification of [26]). A rather complete set of two- loop results for 
the equilibrium dynamics has been only recently reached (see [75] for a brief review). 

In the following we introduce the dynamic universality classes that we will consider 
in this review. 



3.4-1- Purely dissipative relaxation: Model A. The purely dissipative dynamics of a 
iV-component field (fii (i — 1...N) can be specified in terms of the stochastic Langevin 
equation 

fito(x, t) = -^t^4t + ^ x > *) > ( 33 ) 
difi(K,t) 

where H[(p] is the static reduced Hamiltonian of the model, fl is a kinetic coefficient and 
Ci(x, t) a zero-mean Gaussian white noise with correlations 

(Ci(x,i)Ci(x',0) = 2n6(x-x!)6(t-lf)6 ij , (34) 

where T is the bath temperature. Equation (33) is a special case of the more general 
one introduced in section 3.3, with V — fi, Af — Q and it fulfills Einstein's relation 
V = AT. 

The critical dynamics of some anisotropic magnets and alloys [26] are described 
by equations (33) and (34) with H given by equation (22) with N = 1, i.e., by the 
effective Hamiltonian of the Ising universality class. Model A describes also the critical 
dynamics of kinetic spin models on the lattice and spin-flip sampling, i.e., in which the 
elementary step amounts to an arbitrary change in the orientation of the spin in a given 
site, performed with proper rates (e.g., the well-known Glauber dynamics [76]). In the 
case of A-component spin models on a regular lattice and 0(iV)-symmetric short-range 
interactions, the critical dynamics is described by equation (33) where H is given by 
equation (22). 

The dynamical properties of Model A may be worked out by representing the 
Langevin equation (33) as a dynamical functional, following the method outlined in 
the previous section. The resulting action is 



dt J d d x 



<fd t (p + VLip — <pQ<p 

dip 



(35) 



To specify completely the dynamics one has to provide the initial condition for the 
field (p(x.,t): (p(x.,to) = y?o( x )- More generally one can assign a probability distribution 
function for the initial condition, in the form exp{— Tiofv^o]} and then average over the 
initial field y?o( x )- If the system is already in thermal equilibrium at time to, then 
T'My^o] = ^[v^o] and one can equivalently extend the time integration in S[ip, </?], from 
— oo to oo [24]. This is possible since, assuming ergodicity, independently of the initial 
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condition in the far past, the same stationary order parameter distribution is reached 
at time to- The resulting theory is translational invariant both in space and time, and 
given that the Einstein's relation is fulfilled, equal-time correlation functions can be 
computed directly using the functional distribution e~ n ^. 



3.4-2. Conserved order parameter: Model B. In some cases Model A is not suited to 
describe the dynamics of physical systems. For instance, when the order parameter is 
related to the density of particles in a fluid, as it is the case when studying the liquid-gas 
critical point, one expects the continuity equation to be obeyed. The simplest model in 
which such a local conservation law is implemented is known as Model B, according to 
the classification of [26]. In particular, the dynamics of the scalar order parameter ip is 
given by 

3^(x,*) + V-J(x,i) = 0, (36) 

where J(x, t) is a fluctuating driving force. Thinking of (p as a particle density one 
expects the deterministic driving force to be related to the gradient of some sort of local 
chemical potential. Indeed it is generally assumed that 

J(x, t) = -aV x ^M- + J c (x, t), (37) 

where a is a kinetic coefficient and a Gaussian zero-mean random current with 
correlations 

(J c ,,(x,t)J (x',t')> = 2a5(x-x')5(t-t')5 ij • (38) 

In the example given above, a is related to the mobility whereas /2(x, t) = 5H[(p]/8(p(-x, t) 
represents a sort of chemical potential whose inhomogeneities drive the diffusion of the 
particles. 

It is not difficult to recast the dynamic model specified by the previous equations 
in the form of the Langevin equation (33), with the noise correlation (34) [where 
( = -V • J c ], N = 1 and Q -> -aVl [26, 30]. 

Model B describes the critical dynamic properties of some uniaxial ferromagnets 
[26] which belong to the Ising universality class, i.e., the effective Hamiltonian 7i in 
equation (37) is given by equation (22) with N — 1. Moreover Model B describes the 
critical behaviour of lattice spin models with spin-exchange sampling (also known as 
Kawasaki dynamics [77]), characterized by an elementary step which amounts to an 
exchange between the spins in two neighbouring sites. 

The dynamical functional associated with Model B is given by 



S[<p, <p}= dt d d x 



2<^M , ~ V72~ 



The initial condition can be accounted for as described in the case of Model A. 



(39) 
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3.4-3. Coupling of a conserved density to the non-conserved order parameter: Model 
C. The models of dynamics considered in the previous subsections involve only the 
time evolution of the order parameter tp, subjected to the thermal fluctuations due to 
the coupling to the thermal bath. These models, although quite rich and complex in 
their phenomenology, are usually too oversimplified to be able to capture the relevant 
features of real systems, in particular those with a vector order parameter. In many 
cases the order parameter is not the only relevant slow variable that has to be taken into 
account. For instance, in the case of a one-component fluid close to its critical point, 
the conserved order parameter interacts with three slow hydrodynamic modes [31]. 

The simplest model with two interacting fields, called Model C, consists of a 
non-conserved iV-component vector order parameter (p(x.,t) coupled to a non-critical 
conserved density e(x, t). 

The dynamics is specified by the following coupled stochastic Langevin equations 

MM)= -n^^ + c„(x,f), (40) 

d t e(x,t) =^^! + Ce(x,0 , (41) 

where Q and p are the kinetic coefficients, 7i[(p, e] the reduced Hamiltonian of the system 
and C^( x , t), Ce( x )0 zero-mean Gaussian noises with 

(C v (x,*)C v ,(x',*')) = 2fi<$(x-x')<S(t-t')> (42) 
(C £ (x,t)C £ (x',t')> = -2pfiV^(x-x')5(t-0- (43) 

Note that Einstein's relation is fulfilled by this system of coupled Langevin equations, 
therefore the equilibrium distribution function of the fields e and <p is given by 
exp{— H.[<p, £]}, suitably normalized. 

In the case of a system with short-range 0(iV)-symmetric interactions, the 
Hamiltonian H[p,e] in equation (40) is given by 



H[<p,e] =H[p] + J d d x 



1 2 1 2 



(44) 



where TC[p] is the 0(iV)-symmetric Landau-Ginzburg Hamiltonian reported in equation 
(22) and 7 the coupling constant between y?(x, t) and e(x, t). In particular this 
coupling does not change the critical static properties of the field p as it can be seen by 
computing (through a Gaussian integration) the effective Hamiltonian H e ff[¥>] defined 
by exp{—7i c fr [(/?]} = f[de] exp{— H[p, e]}. Indeed it turns out that 7Y c ff[y?] = 'H[ { p] (see 
equation (22)) where go i— > g' = go — 127q- On the same footing one can show that 
e- field equilibrium correlation functions are related to y? 2 -field correlation functions [30]. 
For this reason e is also referred to as energy density, given that p> 2 is conjugate to the 
temperature (oc r ) in equation (22) and therefore the derivatives of the free-energy with 
respect to the temperature give rise to y? 2 (e)-correlation functions. 

Some lattice models belonging to this universality class are discussed in [78]. 
Among real systems whose dynamic critical properties are described by Model C with 
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Hamiltonian (44) and N = 1 we mention the case of structural (displacive) transitions 
in crystalline materials and some phase transition in uniaxial antiferromagnets [26]. 

Dynamical correlation functions, generated by the Langevin equations (40) and 
averaged over the noises ( v and ( £ , may be obtained by means of the field-theoretical 
action 



S= / d d xdt 



0dt<p + £l<p SH[ ? ,£] - ^ + ed t e - P niVl 6H[ f ,£] + epQV 2 J 



*p — — - ^ Se — j' (45) 

where two response fields </3(x, t) and e(x, t) associated with (p(x.,t) and e(x, t) have 
been introduced. It is easy to read from equation (45) and (44) the interaction vertices, 
given by — ilg (f>(p 3 /3\ (as in the case of Model A), — Q^f e(p(p and ptl'fo (p 2 V 2 e/2. 

To take into account the effect of the initial condition one has also to average over 
the possible initial configurations of both the order parameter y?o( x ) = ty?(x, t = 0) and 
the conserved density £o( x ) = e(x,t = 0) with a probability distribution e - H o{<po,eo] gj ven 
by [81] 

H [po] = Jd d x {f h(x)- M (x)f + ^[ £o (x)-<x)f}. (46) 

This specifies an initial state w(x) for y?(x, t) and v(x) for e(x, t) with correlations 
proportional to Tq 1 and Cq, respectively. 

The model with a conserved order parameter coupled to a conserved density is 
called Model D [26]. 

3.4-4- Other dynamic universality classes. The models of critical dynamics introduced 
so far are in many cases still too simple to describe the critical dynamic properties of a 
variety of actual materials. 

For instance let us briefly discuss the case of isotropic ferromagnets, microscopically 
described by the so-called Heisenberg model for the three-component spin variables 
(see, e.g., [28]). In the absence of symmetry-breaking fields the model has a 0(3) 
symmetry. At a mesoscopic level the system is described by a three-component 
order parameter S(x, t), representing the local magnetization, which turns out to be a 
dynamically conserved quantity. In view of the underlying 0(3) symmetry, the universal 
aspects of the equilibrium critical properties of the model are correctly captured by the 
Hamiltonian (22) with N = 3 [67] (ipi h- > Si). The mesoscopic dynamics of order 
parameter S(x, t) has to account both for the conservation law mentioned above and for 
the expected Larmor's precession of S(x, t) in the local magnetic field 5Tl[S]/5S(x.,t) 
generated by neighbouring spins. The Langevin equation accounting for these features 
and the associated universality class are known as Model J [26]: 

a t S(x,i) = - 9l S(x,i) x Jj^L + aVl-^K. + C ( Xj () , m 

where gi is the coupling constant of the precession term, a the kinetic coefficient 
(analogous to that of Model B) and £ a Gaussian random current with correlations 
given by equation (34) where Q \— > -crV^. Apart from the precession term (which is 



Ageing Properties of Critical Systems 



23 



an example of the general class of reversible mode- couplings, see [28, 30] for details), 
the dynamics of each single component Si is given by a Model B and satisfies Einstein's 
relation. The precession term, on the other hand, does not influence the equilibrium 
distribution which is still given by exp{— 7Y[S]}, properly normalized [28, 30]. We do 
not discuss the properties characterizing this universality class but we mention only that 
they agrees with what has been observed in experiments [26]. 

The critical properties of an isotropic antiferromagnet close to the Neel temperature 
are instead described by a non-conserved three-component order parameter m s (x, t) (the 
local staggered magnetization) which is dynamically coupled, via precession-like terms, 
to a three-component conserved density representing the local magnetization m(x, t). 
The resulting dynamic universality class is known as Model G and its predictions 
agree with the experimental results on actual magnetic systems. We remark that the 
equilibrium properties of the model are given by the Hamiltonian (22) (with N = 3) for 
the order parameter m s [67]. Accordingly, models G and J provide a non-trivial example 
of two different dynamic universality classes of real magnetic systems that reduce to the 
same static universality class. Other examples are provided by the dynamic universality 
classes describing a planar magnet and the superfluid 4 He, known as Model E and F, 
respectively [26], which belong to the same static universality class as the Hamiltonian 
(22) with two-component order parameter (i.e., 0(2) symmetry). 

Since we do not consider these models in what follows, we do not discuss any further 
details of the associated dynamic universality classes and refer the interested reader to 
the available review [26] and textbooks [31, 28] on the subject. 

3. 5. Renormalization and scaling of a purely dissipative critical system 

As we discussed in subsection 2.1 the quantities one is interested in when studying 
dynamical processes are the response and correlation functions. In the following we 
will generically refer to them as correlation functions. In subsection 3.3 we described 
the method that allows the computation of such functions as averages of the form (29) 
where the specific expression of the functional S to depends on the considered dynamics. 
Field-theoretical methods can be applied to compute such averages. As usual, when 
going beyond the Gaussian approximation (i.e., considering the anharmonic terms in 
the expansion of S to ) one faces the problem of ultraviolet (UV) divergences that can be 
regularized by standard methods and removed by means of a standard renormalization 
procedure. The existence of a well-defined renormalized theory allows the derivation 
of RG equations for the correlation functions. In particular their solutions highlight 
the scaling properties of these functions in the large-distance long-time limit [the so- 
called infrared (IR) behaviour] which can be computed by means of the renormalized 
perturbation theory. 

To illustrate the general procedure, we consider as a specific example the Model A 
dynamics of a ip 4 field theory. All the following basic steps can be carried over to more 
complex cases. 
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3.5.1. RG, e expansion and Minimal Subtraction. The RG idea is to construct a 
scale-dependent effective action S^[ip, (p\ (with x > 1) which has connected correlation 
functions G n ^ = ([tp] n [fi] n ) c satisfying 

G$({x, 0) = Z- nt \K)Z^I\K)G n ^{^ xH}) + ..., (48) 

where is computed with the effective action S w whereas G n ^ is obtained with 
the original action S. One can arbitrarily fix S^=i = S. The ellipsis in equation (48) 
stand for functions decreasing faster than any power of x for x — > oo, {x, £} stands 
for the collection of the n + n space-time points, and we allow an anisotropic scaling 
between space and time through the exponent z. The mapping S[ip, Cp\ — > ^[y, is 
called RG transformation and {S^ty?, i is called the trajectory of S[tp, (p\ under 
the RG flow. Various RG transformations differ by the form of Z(x), Z(x), and of 
Sx- In explicit constructions the transformations are generated by integrations over the 
large-momentum fluctuation modes of the fields. 

The coupling constants appearing in S^ip, (p\ are explicit function of x. If for 
x — > oo, the action S^[(p, 0\ has a limit S*[<p, </?], called the fixed-point action, the 
correlation functions behave like 

G niJi ({xx,x^})~Z:/ 2 (x)ZlV 2 (x)G; ift ({x,t}) for x ^ oo , (49) 

where * is used to denote fixed-point quantities and G* n fl is computed with the fixed- 
point action S*. It is possible to show that under some reasonable assumptions [30] 
Z*(x) = x -2 ^ and >2*(x) = x~ 2c ^. Accordingly, correlation functions have a scaling 
behaviour at large distances and long times: 

G n ,n({xx, H*t}) ~ x-^-"^G; ift ({x, t}) for x -> oo , (50) 

where the exponents d v , (the so-called scaling dimensions of the fields), and z are 
properties of the fixed point, i.e., they depend only on S*[ip, Cp\. Correlation functions of 
all those actions that flow under RG transformations into the same fixed-point action S* 
display the same critical scaling behaviour: This provides the basis for the universality 
observed in critical phenomena. Within this framework a given universality class is 
characterized by the values of the exponents d v , dq, and z and by the scaling function 
G* n h . There are several ways and approximations that allow to calculate such quantities 
in physical dimensions d — 2, 3. 

Renormalization-Group equations for the correlation functions can be conveniently 
derived within the field-theoretical approach. They eventually lead to a relation of the 
form (48) for correlation functions. The first step is to regularize the theory in such a 
way to render finite those Feynman diagrams of the perturbative expansion that would 
be otherwise UV-divergent. To this end we will assume dimensional regularization [30]. 
The original UV divergences of the four- dimensional theory appear in the correlation 
functions as poles when e — > 0, i.e., when the regularization is formally removed. Let 
us consider the Model A dynamics with Hamiltonian given by equation (22). In this 
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specific case one introduces the renormalized parameters as: 

° ° -1 ° 4— d 

Q = r = r = Z Z r r, and g = g = Z g n g. 

Here the parameters with the superscript o are those originally appearing in equation 
(35) {bare parameters) whereas the others represent the renormalized ones. Note that 
the renormalized coupling constant g is defined in such a way to be dimensionless. \i is an 
arbitrary scale with the dimension of a mass that in general enters in the renormalization 
of the parameters that are not dimensionless. The Z factors are called renormalization 
constants and can be computed in many different ways. In the following we will use 
the so-called Minimal Subtraction (MS) scheme [30] which allows one to study directly 
the critical theory, hence we fix r = in what follows. Within MS the renormalization 
constants take the form Z — 1 + ^2oi n (g)/e n , where the coefficient a n (g) = 0(g n ) are 
computed order by order in perturbation theory. It is possible to show that the FDT 
yields the relation Z w = (Z/Z) 1 ' 2 [30]. O nee the correlation functions of the fields are 
expressed in terms of the renormalized parameters only, they turn out to be regular 
functions for e — > 0. 

To renormalize G n ^ one uses the Z-factors in equations (51) obtaining 

Gn,n — Z n / 2 Z n / 2 G nj n , (52) 

where the r.h.s. is expressed in terms of renormalized quantities. The RG equations may 
be derived by exploiting the fact that the bare correlation functions are independent of 
the mass scale \i introduced to define the renormalized theory 

o 

^ G n ,h = , (53) 

where the notation | reminds that the derivative is taken with fixed bare parameters. 
As we show below, this equation defines a RG flow in the parameter space {g, /i, Q}. 
In terms of renormalized quantities, equation (53) reads 



d n„ n d d 

,^ + 2 7+ 2 7 + 7 ^ + ^. 



G nifi ({x,i};<7,A*,ft) = 0, (54) 
where we the so-called Wilson's functions are defined by 

l{g) = ^lnZ| , 7(0) = //^lnZlo, 7 W (0) = /i<9 M lnfi| , (55) 
and 

/%)=M^| , (56) 

which are regular for e = [30]. 

To solve RG equation (54), the method of characteristic is usually employed. One 
introduces a continuous real scale parameter /, in terms of which the flowing mass scale 
/i reads = \il and \i plays the role of an initial condition for I — 1. The IR behaviour 
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is obtained for / — > 0. Then one defines effective /-dependent parameters by solving the 
ordinary differential equations 

z^ = n(0 7w (0, 1^ = 0(0, etc., (57) 

with f2(Z = f) = fi, g(/ — 1) — g etc., where the RG functions on the r.h.s. depend on 
I through g(l), e.g., j(l) = / ~f(g(l)). The RG equation for the renormalized correlation 
functions becomes an ordinary differential equation 



j7W + j7(0 + ' S 



G„,»({x,t}; 9 (/),M0.n(0) = C (58) 



whose solution is 

G n>fi ({x, 0; MO, m) = Gn,n({x, 0; <?, exp + W)] } • (59) 

Physically, varying the parameter /, one changes the scale at which the model is 
explored. The infrared regime is reached for I — > 0. The flowing parameters 0,(1) and 
g(l) can be interpreted as effective values of the original ones on different length/time 
scales set by fil. They evolve under scale transformation \x — > fil according to the flow 
equations, that are determined by the required renormalizations through the functions 
7, 7, 7w and (5 (see equations (55) and (56)). 

Under a scale transformation, the dimensionless renormalized coupling g changes 
according to equation (56). Therefore g(l) increases for I — > when the f3 function is 
negative and decreases in the opposite case. Fixed points for the flow of g correspond 
to zeroes g of the (3 function, i.e., (3(g) = 0. Those where (3'(g) < are not stable in 
the limit I — > (IR behaviour) we are interested in: The effective coupling moves away 
from them. Conversely, when f3'(g) > 0, the zero is stable. Indeed linearizing the flow 
equation for g close to g one finds that g(l) — g = (g — g)l^'^ for I <C 1. Let us assume 
we find such an IR-stable fixed point for g = g*. We denote the fixed-point values of 
the RG functions as r\ = j(g*), rj = 7 (g*), r) u = 7w(<7*) = (v — ? ?)/2 (this last equality 
follows from the FDT), and z = 2 + i] w . Accordingly, Q(l — > 0) = Ql z ~ 2 . 

The leading IR behaviour is obtained from equation (59) 

G n , fi ({x, 0; g, -i^o r i/2+ ^ /2 G n ^({x, 0; 9(i),^(i), m) ■ (60) 

From the dimensional analysis one finds that 

G B>fi ({x, 0; 9, PL, SI) = i-^-^'^^G^/i, t/r}; g, tf, Or/f) , (61) 

where £ and r are arbitrary units of measure of length and time, respectively. Applying 
this equation to the r.h.s. of equation (60) with £ = (^l)' 1 and r = VL^ 1 i 2 l 2 ~ z : , one 
finally arrives at the leading IR scaling behaviour of the correlation functions 

G n ,n({x, 0; 9, /i, n) = Z^ fi )/i B ^ 2 )/^*- 2 )/ 2 G B , fi ({Z0ix), f)}; <?*, 1, 1) , (62) 

where 

5(n, h) = n ~*~ - + n — — . (63) 
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In the following, to simplify the notation, we will set the constants \i and Q to 1. 

Let us consider in more detail the scaling forms of the two-point critical correlation 
function C q (t, s) and response function i? q (t, s), with t > s. From equation (62) and 
equation (63), using TTI, one has 

C q (t,s) = Jd d xe^G 2fi = l^Mr^nt - s)) , 

R^(t,s) = n Jd d xe i ^G 1 , 1 =l r ' +z ^f R (r 1 q,l z (t-s)) . (64) 
If we choose I = (t — s) -1 ^, with i - s » 1, these scaling forms may be rewritten as 

(? = |q|) 

C^s) = (t- s)^ z F^(q(t- s) 1 '*) , 

R^t, s) = (t- sf'^F^t - sf/ z ) , (65) 

where F^(x) and F^(x) are regular functions for small argument. The FDT provides 
a differential relation between F^(x) and F^(x). 

By means of RG, one reduces the problem of calculating universal quantities to the 
evaluation of RG functions at the IR-stable fixed point. In order to find such a fixed 
point, a useful tool is the e expansion. It is based on the observation that d = 4 is a 
special dimension for the model we are considering, in fact using the definition of the f3 
function we get 

M = M = "(4 - d) {^f^j 1 • (66) 

Being Z g = 1 + 0(g), the fixed point g* = is IR stable for d > 4 and unstable in 
the opposite case. Thus, for d > 4 the critical behaviour is governed by the fixed-point 
g* — (Gaussian fixed point), corresponding to the Gaussian model. For d < 4 this is 
no longer true. However, adiabatically decreasing the dimension from four dimensions, 
the new IR stable fixed point is expected to be close to the Gaussian one, i.e., g* = 0(e), 
where e = 4 — d. As originally pointed out in a seminal paper by Wilson and Fisher [79], 
one can perform a double expansion of RG functions in terms of g and e and find the 
zero(es) of the (3 function as series in e. Finally the critical exponents and other universal 
quantities are obtained expanding in e the corresponding RG functions evaluated at the 
stable fixed point. 

This procedure allows one to obtain the critical quantities as series in e. The 
analytical continuation of such series at physical dimensions d = 3, 2 (i.e., e = 1,2) 
may seem not so straightforward. However, the validity of this continuation to 
e = 1,2 is corroborated by the very good agreement of the e-expansion estimates 
(obtained applying resummation techniques based on Borel resummation [30]) with 
other theoretical and experimental values [67, 30]. 

3.5.2. N on- equilibrium renormalization and scaling. Let us consider the case with a 
non-equilibrium initial condition. As explained in section 3.4, it can be studied by 
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averaging the correlation functions with the probability distribution function of the 
initial condition, specified by exp{— 7i } [24]. In turn this implies that expectation 
values in equation (29) have to be computed with a total dynamical functional given 
by S[tp, 0\ + Tio (see equation (35)). As long as Tio has the same form as the static 
Hamiltonian 7i, but with different bare couplings, the usual renormalizations are enough 
to render the theory finite [24, 25]. This is no longer true if, for instance, the initial 
state is an uncritical one like a high-temperature state with short-range correlations and 
a small initially prepared magnetization (</? (x)) = a(x), i.e. 

(Mx) - a(x)][y. (x') - a(x')]> = r^x - x') . (67) 

The corresponding 7io[y?o] is Gaussian 

Wobo] = |d^Mx)-a(x)] 2 . (68) 

The canonical mass dimension of r is [r ] can = 2, thus the possible fixed points of r 
can be ±oo and 0. Since r = and — oo yield non-normalizable distributions, the fixed 
point of physical interest is r = oo [24]. This corresponds in the language of surface 
critical behaviour to the ordinary transition [80]. Corrections due to finite value of t 
will be irrelevant in the RG sense, thus one can set Tq" 1 = from the very beginning of 
the calculation, as long as the leading contribution does not vanish with Tq 1 = 0. 

Following standard methods, the response and correlation functions may be 
obtained by a perturbative expansion of the functional weight exp{— S[ip, <p\ — 7i [y?o]} 
in terms of the coupling constant go- From the technical point of view, the breaking 
of TTI does not allow the factorization of connected correlation functions in terms of 
one-particle irreducible ones, as usually done when TTI is not broken. As a consequence 
all the perturbative computations must be done in terms of connected functions only 
[80, 24]. 

The addition of 7io[y?o] gives rise to new divergences in perturbation theory 
whenever s approaches the "time surface" located at to = 0. On the other hand, 
the bulk renormalization functions defined by equations (51) are not changed by the 
presence of 7io[y?o]> if MS is employed to renormalize the theory. It is possible to remove 
these new singularities by counterterms "located" at that time surface, as it has been 
shown in [80] for the case of surface critical phenomena. Moreover, being the canonical 
mass dimension of time variables —2, the degree of divergence of a generic correlation 
function decreases by 2 for each vanishing time argument. As a consequence the new 
renormalizations are required only in the case of two-point functions. A detailed analysis 
shows that once equilibrium theory has been renormalized according to equation (51), 
only one new renormalization constant is required to render finite both the correlation 
and the response functions. In the case of non-equilibrium Model A dynamics the new 
renormalization is 

^o(x) ^ (Z Z) 1 /^o(x), (69) 
and Zq is the new renormalization constant. 
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The scaling properties of connected correlation functions 

G% = (bn^oP) (70) 

may be exploited by using RG equations as in equilibrium. To renormalize G™^ one 
uses all the Z-factors previously introduced (equations (51) and (69)), obtaining 

o no _ _ 

G% ^ G n , h = Z^Z^{ZZ Q f^G% , (71) 
where again the r.h.s. is expressed in terms of renormalized quantities. The RG 

= 0, in analogy with equation (53). To solve 



o «-0 



equations are derived from /i<9 u G n n 
it, a new Wilson's function 70 = /i<9 u 



nZ |o has to be introduced. The leading scaling 



behaviour of the correlation functions (70) is determined as in equilibrium 

G?„({x, t}) = 1 5 ^^G%{{1^ Ft}) , (72) 

where now 

x , ~ ~ , d-2 + r, „ d + 2 + ft . d + 2 + f, + r, 
b(n,n,n )=n - hn - h n - , (73) 

with the novel exponent i] = 7o(g*), in terms of which the so-called initial-slip 
exponent 6 = —r) /(2z) is defined [24]. The scaling dimensions of the surface field 

is [<A)]scal = (d + 2 + fj + 7/o)/2. 

According to the same procedure, the analysis has been done for various models. 
Model C dynamics (see sections 3.4.3 and 7) has been studied in [81], whereas Models E 
and G [see section (3.4.4)] have been studied in [82]. Model A dynamics of a tricritical 
point is investigated in [83] . In [84, 85] the weakly dilute Ising Model has been considered 
in the case of uncorrelated impurities. The case of extended random defects has been 
also addressed [86]. Other interesting universality classes have been studied [87]. 

Let us consider in more details the scaling forms of the two-point critical correlation 
function C q (t, s) and response function i? q (t, s), with t > s > (i.e., with t and s in the 
"bulk"). From equation (72) and equation (73) one has (we set Q — 1) 

C q (t,s) =G° = r- 2 C z -i q (Z*t,Z*s) , 

R^t, s) = G° ! = Z"+*- 2 i2,-i q (Z**, l z s) . (74) 
Considering I = (t — s)^ 1 ^, these scaling forms may be rewritten as 
C q (f , s)=(t- sf-^ z F c (q(t - sfl\ s/t) , 

i? q (t, s) = (t- sf-^ z F R (q(t - s) l l\ s/t) , (75) 

where, as a fundamental difference to equilibrium ones, they depend on the ratio s/t. 

We observe that the scaling functions Tc and Tr just introduced are not expected 
to be regular when s approaches the time surface, i.e., for s — > 0. To determine the 
functional form of the correlation functions when their arguments approach exceptional 
points, the short-distance expansion (SDE) must be used (see [30, 57] for general 
reference and [80] for applications to surface critical phenomena). The starting point 
is a formal expansion of the fields y?(x, s) and y5(x, s) around s = 0. First of all one 
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notes that, when inserted in correlation functions with bulk fields, the following relations 
hold [24, 25] 

</?(x,0) = <p (x) = and <9 s y?(x, s)\ 8=Q = <p (x) = 2fi<^ (x) . (76) 

As a consequence, for small s, one can formally expand the fields as 

y?(x, s) ~ 0(s)(^ o (x) + h.o.c.f. , (77) 

<p(x, s) ~ 0(s)(^ o (x) + h.o.c.f. , (78) 

where h.o.c.f stands for higher-order composite fields which are neglected if one is 
interested only in the leading contributions in the limit s/t — > 0. Inserting the 
relations (77) and (78) into the correlation functions and taking into account the scaling 
behaviour equation (72), one deduces that, at criticality [24, 25], 

0(s) = acs 1 ' 9 , 

4>(s) = a R s- e , (79) 

(i.e., [0] sca i = [y?] sca i - fiscal = -z(l - 0) and [0] scal = [(^] scal - fiscal = zO, with 
[time] sca i = —z) where ac and a R are two non- vanishing constants. Accordingly for 
small s one has 

C q (t,s) = (f)(s)(ip(q,t)<p (-q)) , 

R^(t,s) =^)Mq,^ (-q)) • (80) 

From equation (72) the scaling form of (ip(q,t)<fo(—q)) = G\ ({q,t};r) can be 
determined and at criticality 

G li0 ({q, 0; o) = r^-'-^ottr 1 ^ o) . (si) 

Taking into account the previous three relations, the following conclusion may be drawn 
(?=|q|) 

C q (t,a) =a c t^'(t/s) e - 1 fc(qt 1/z ) , 

R^t, s) = a R t^^l z {t/s) d f R {qt l l z ) . (82) 
Comparing these forms with equation (75) one concludes that for y — > 
Fc(x,y) ~ a c y' e+1 fc(x) , 

F R (x,y)~a R y- e f R (x) . (83) 

It is possible to rewrite equation (75) in terms of scaling functions F c (x, y) and F R (x, y) 
with a regular behaviour for y — > 0, i.e., 

C q (f, s) = (f - s^^'it/s^Fdqit - s) l l-\ s/t) , (84) 
i? q (i, s) = (t - 8)^^/8)° F R (q(t - s) 1 ^, s/0 . (85) 

These results refer to systems in the infinite- volume limit. Finite-size effects on the 
non-equilibrium dynamics have been studied in [88]. 

For q = the scaling forms (84) and (85) may be written as 

R n=0 (t, s) — A R (t — s) a (t/s) e F R (s/t) , (86) 
C q=0 (t, s) = A c s(t - s) a (t/s) e F c (s/t) , (87) 
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where a = (2 — rj — z)/z. We singled out explicitly the non-universal amplitudes A R c 
by fixing F R; c(0) = 1- With this normalization, F R; c are universal scaling functions. 
From the previous scaling forms one deduces that 

dsC^t, s) = A dc (t - s) a (t/s) 9 F dc (s/t) , (88) 

where the non-universal amplitude Age has been defined in such a way that Fgc(0) = 1. 
The relation between Fq C (x) and Fq(x) can be easily worked out. The amplitudes are 
related by A dc = A c (l - 9). 

Using equation (86) and (88) one finds that 

iVo(M) A R F R {s/t) 
q=0i ' } ~ dsC^t, s) ~ A dc F dc (s/t) m 
is a universal scaling function, being the ratio of two quantities (R q=0 (t,s) and 
(9 s C q=0 (t, s)) that have the same scaling dimensions. Furthermore it is a function only 
of s/t and not of s and t separately. 

From equation (89) and using the results of section 2.4, one easily finds that 

X°° = lim lim - R q=°( t ; s ) = i im XJs/t) = 4 s - = - A AR tt, , (90) 

is a universal, dimensionless amplitude-ratio in the sense of [69]. 

Fourier transforming equations (84) and (85) we obtain the scaling forms of the 
autoresponse and autocorrelation functions 

ik=o(t, s) — A R (t — s) a - d l z (t/s) d F R (s/t) , (91) 

Cx=o(*, s) = Ac s(t - s) a - d l z {t/s) d F c {s/t) , (92) 

with A R , Ac, ^Fr{.u) and J 7 c(y) that are (in analogy with their counterpart in the real 
space) non-universal and universal, respectively. 

Note that even X x=0 (t, s) is a function of s/t only. This is an important difference 
compared to mean-field glassy models, where X x=0 can be written as a function of 
C X=0 (M) (see, e.g., [7]). 

3.6. The effect of the initial condition and the irrelevance of Tq 1 . 

In the previous section, by means of dimensional analysis, we concluded that Tq 1 is 
an irrelevant variable and so it can be neglected as long as one is interested in the 
leading long-time behaviour. However, this is not possible when considering quantities 
that vanish for t^ 1 = 0. This is actually the case of the correlation function when 
the smaller time s goes to zero, as it is clear from equations (87) and (92). These 
forms have been obtained assuming large s, but it is simple to show that due to the 
Dirichlet boundary condition in time (implicit in the assumption t^ 1 = 0) the correlation 
functions with the insertion of one <£>(x, 0) vanish identically, see equation (76). In these 
cases the contribution depending on t^ 1 is the leading one and so the scaling forms can 
not be calculated setting r^ 1 = 0. It was shown that the insertion of a boundary field 
is equivalent to the insertion of r^Vo [25]: 

¥»(x,0) = r o -Vo(x), (93) 
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when inserted in a correlation function with bulk fields. Therefore, using equations (72) 
and (73), the scaling of the autocorrelation function for s = and large t turns out to 
be (neglecting higher-order corrections due to initial correlations, i.e., to finite To): 

A(t) = C x= oM) = r - 1 ( V ,(x,t)^ (x)) ~ To- 1 *-* 1 ' ' 1 '/* = r -Y'-^, (94) 

where we introduced the so-called magnetization initial-slip exponent, given by 6' = 
[d — 5(1, 0, l)]/z [see equation (73)]. For Model A dynamics it can be written as 

0f = e + a = 9+ (2-r,-z)/z. (95) 

Equations (94) and (92) imply that for t > s, C x=0 (t,0) ~ r ~V d/2:+(? ' whereas 
C x=0 (t, s) ~ s i-d-f--d/z+e' ^ ^ leading order in t^ 1 . In both the cases the contributions 
decrease as t~ d / z+e ' , allowing for the identification of the autocorrelation exponent, 
defined as [46] 

C x=0 (t, s) ~ r x/z for t > s , (96) 

so that 

X = d-9'z = d-9z-2 + ri + z. (97) 

Let us point out that the large-t behaviour of C x=0 (£,0) and C x=0 (t, s) can, a priori, 
be different, giving rise to a temporal cross-over for finite s, as in the case of Model B 
dynamics considered in section 6. 

For sake of completeness we comment that also autoresponse exponent is usually 
defined in the literature [89], through 

i?x=o(f , s ) ~ t~ Xa/z for * > s ■ (98) 

For Model A dynamics \ R equals A. This is not the case in the more general case, see 
section 6. 



3. 7. Comments 

Some remarks are now in order concerning the particular choices of the noise and of the 
initial conditions. 

If in equation (33) we fix C( x ?^) = 0? we obtain the deterministic evolution of 
an order parameter under a diffusion-like equation. In turn, the associated dynamical 
action has no quadratic term in (p. Physically this can be used to describe the evolution 
of a system quenched from a fully disordered initial state to a critical point, at which the 
thermal noise is negligible. This should be the case of systems with T c = 0. However, 
it is not clear under which conditions the microscopic models defined via a master 
equation (and hence stochastic) are in the same universality classes as the corresponding 
deterministic mesoscopic equations. Examples of different resulting universality classes 
have been provided [90]. The RG study and e expansion of this deterministic model 
has been worked out in [66]. For our discussion, a relevant finding is that response and 
correlation functions scale in time with A 7^ A#, resulting in a trivial FDR. For this 
reason we will not further consider the deterministic evolution. 
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We have also assumed the noise to be Gaussian and white. A non-constant spectral 
density of the noise may also occur (in this case the noise is called "coloured"). To our 
knowledge, the non-equilibrium behaviour under the effects of a coloured noise has been 
considered only for the random walk [91], corresponding to the Gaussian approximation 
of the q = mode in equation (33). 

So far, we considered only the case of short-range correlations in the initial state. 
Long-range ones (for instance decaying as ^ | x l _p f° r large x) may change the dynamic 
non-equilibrium universality class and in particular the initial-slip exponent 9 and X°°. 
Decreasing the exponent p of the long-range correlations, the canonical dimension of 
To, that is negative for large enough p, may become positive. In this case 9) turns to 
a relevant parameter in the RG sense and one expects that the stable fixed point and 
consequently the critical properties get modified. Exact calculations for the spherical 
model [89] and for the one-dimensional Ising model [92] substantiate this picture. 

Finally a scaling correlation function of the form of equation (92) was also found 
in models of self-organized criticality [93]. 



3.8. Scaling of high- order observables 

In this section we briefly discuss the scaling behaviour of observables that are powers of 
the order parameter at the same space-time point (composite operators in FT language). 
More details can be found in reference [23]. 

Let us consider an observable O with ho as a conjugate field (e.g., O is the energy 
density and h the temperature), that couples to H according to H i— > H — (3 h 0. As 
a consequence, the dynamical functional S changes according to S \— > S — f3 hoO, where 
the associated operator O is given by 



O 



y did ,^ (x , ()I ,_j|_. (99) 

The linear response of an observable A to a variation in the field ho conjugated to O 
can be expressed as 



5h 



o 



= (3 {AO), (100) 

h o =0 

where {-)h stands for the average over the dynamics associated with the dynamical 
functional in the presence of ho- 

We are interested in the scaling properties of the correlation functions of the form 
(OO) and (OO). Using RG equations as in the previous sections one obtains [23] 

(0(q, t)0(-q, s)) = (t- sT +1 F c (q z (t - s), s/t) , 

(0(q, t)0(-q, s)) = (t- sy°F R (q z (t - s), s/t) , (101) 

where ao = (2 — rjo — z)/z and rjo is the scaling dimension of the observable O. 

In the presence of a time surface the functions F c {q z {t—s), s/t) and F R (q z (t—s), s/t) 
are not regular in the limit s/t — > 0, as in the case of O — tp of the previous section. To 
extract the leading singularity in the ageing regime a proper short- distance expansion 
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should be considered. This can be easily worked out for 0^ m \t) = y? m (t) (hence 
0^ m \t) ~ &(t) V m - l (t)). In this case the analysis is straightforward and the final scaling 
forms are [23] 

(C (m) (q, t)0 (m) {-q, s)> = s(t - s) a °{t/s)- (m -V +me Fg{q z {t - s), s/t) 
(O im \q,t)d (m \-q,s)) = {t- s) ao (t/s)- {m ~ 1)+me F°(q z (t- s),s/t) (102) 

where the functions F® and F R are regular for s/t — > 0. Furthermore they are also 
universal once we fix the normalization for small arguments. 

Let us focus on the scaling properties of the correlation and response functions of 
the two zero-momentum observables with m = 2: 

N N 

E(t) =J2 ddx v^ 2 (x, t) = J2 ( d <?) ^(q, *)?i(-q> > ( 103 ) 
i=i ^ i=i ^ 

TW0= Jd d x<p i (x,t)<p j (x,t) = J{dq) ¥»i(q,t)^-(-q,t) . (104) 

and T^j are usually referred to as the "energy" and the quadratic tensor operator, 
respectively [23]. The corresponding response operators are 

N 

E{t) = (dq)2 l p i (q,t)cp i (-q,t) 
i=i ^ 

TWO = y (dg) [&(q, t)cp 3 (-q, t) + <^(q, *)&(-q, t)\ . (105) 

In the following we will generically refer to them as O(t). In this case, being the scaling 
dimension of 0(t) well-known, the exponent ao can be written in terms of the usual 
critical exponents [23]. The corresponding scaling forms for the correlation and response 
functions (the latter defined as in equation (100) with A — O) are 

C°(t, s) = (0(t)0(s)) = A° c s(t - sy°(t/sr 1+2d F$(s/t) , (106) 
R°(t,s) = (0(t)0(s)) = A®(t — s) a °(t/s)~ 1+2e F®(s/t) , (107) 

where the non-universal amplitudes Aq R are chosen in such a way that Fq R (0) = 1. 
(As previously done, we have absorbed the factor [3~ l in the definition of R°(t, s).) In 
terms of these quantities we can write the (zero-momentum) FDR as 

*^ s 5& (io8) 

that, as its analogous (89), depends only on the ratio s/t, and it is a universal function. 
In particular its long-time limit 

1 A° 

X% = hm lim X {t, s) = - §rjo , (109) 

is universal (analogously to equation (90)). These considerations hold in general for all 
the FDRs that can be obtained from equations (102). 
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4. Purely dissipative dynamics of an O(N) model (Model A) 

One of the simplest non-trivial models displaying ageing after a quench to the critical 
point is a lattice spin model in d dimensions with O(N) symmetry, evolving according 
to a purely dissipative dynamics (no conservation laws). In the simplest instance its 
Hamiltonian is given on the lattice by 



where s ; is a iV-component spin located at the lattice site i, with s ; 2 = 1. Here and in 
the following the symbol (ij) means that the sum runs on all nearest-neighbour pairs 
of lattice sites. When N — 1 the Hamiltonian (110) describes the Ising model, whereas 
for N = 2 the XY or planar model and for N = 3 the (isotropic) Heisenberg model. A 
purely dissipative dynamics for the lattice model (110) proceeds by elementary moves 
that amount to random changes in the direction of the spin Si (spin-flip sampling). The 
transition rates can be arbitrarily chosen provided that the detailed-balance condition 
is satisfied. For analytical studies the most suited is the Glauber dynamics [76], which 
allows some exact solutions. Given its relative simplicity, this lattice model is the 
most studied and best understood. In the following we review analytical and numerical 
results that have to be compared with the field-theoretical (FT) ones presented in the 
next sections. In fact the coarse-grained continuum dynamics described by the Langevin 
equations (33) is expected to be in the same universality class as lattice models with 
O(N) symmetry, short-range interactions and spin-flip dynamics [26]. All the reviewed 
determinations of the FDRs are summarized in table 1 (Magnetization FDR) and table 
2 (other observable FDRs) together with the FT estimates. 

One- dimensional Ising model: In this case there is no finite temperature phase 
transitions, since d — 1 is the lower critical dimensionality of the Ising model, therefore 
T c = 0. Analytical solutions for the critical response and correlation functions have been 
provided [94, 95, 96, 97, 58]. From them the critical exponents read z = 2 and 9=1/2 
(as first shown in reference [98]). X°° = 1/2 has also been found [94, 95] (actually this 
result was indirectly obtained several years before in [99]). Correlation and response 
functions of more complex observables have been studied [37, 58, 97] and in particular 
the FDR for the energy has been found to be X% = at T = 0. The fact that X°° ^ X% 
is interpreted in [58] as an interplay between criticality and coarsening, a peculiarity of 
those models with T c = 0. For many other observables, instead, the equality Xq = X°° 
between the critical FDRs holds [37, 58, 97]. In this context it must be recalled (as 
stated in previous sections) that the correspondence between the Glauber-Ising model 
and the Landau- Ginzburg dynamics is not obvious and they probably belong to two 
different universality classes (in fact they display different non-equilibrium exponent 9 



d-dimensional spherical and 0(oo) models: The static critical behaviour of the 
O(N) model in the limit iV — > oo is equivalent to that of the spherical model, defined 




(110) 



[90]). 
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by the Hamiltonian 

<ij> 

where Si are real numbers subjected to the constraint J2i s i = L d (L being the linear 
dimension of the d- dimensional lattice, assumed for simplicity to be hypercubic). Indeed 
Stanley proved [100] that the free energies of the two models are exactly the same. 
From this equality it follows that critical exponents, universal scaling functions etc. are 
equal. The same equivalence also holds for equilibrium critical dynamics defined in the 
spherical model by means of a Langevin equation properly modified in order to prescribe 
a dynamics that is compatible with the spherical constraint. As far as we are aware, a 
proof of this equivalence has not yet been carried over to the non-equilibrium critical 
dynamics we are interested in. Let us recall that the free energies of the models in the 
presence of a spatial boundary (whose corresponding field theory looks very similar to 
the non-equilibrium one that we are considering) are not equal [101]. 

The spherical model has attracted a lot of attention, since its essentially Gaussian 
Hamiltonian makes it exactly solvable, even though the resulting critical behaviour 
is not mean-field-like in 2 < d < 4, because of the spherical constraint. Even in 
the O(oo) model several quantities are exactly calculable because the fourth order 
interaction term in the Hamiltonian (22) can be self-consistently decoupled according 
to g/N((p 2 ) 2 — > gC^ =0 (t)ifi 2 and the resulting Langevin equation (33) is linear and thus 
solvable [24]. The dynamical critical exponents for both O(oo) [24] and spherical models 
[17] are z = 2 and 6 = 1 — d/A, for 2 < d < A. The spherical FDR of the magnetization is 

= l — 2/d [17]. This result has been found in agreement with the direct calculation 
for the O(oo) model [23]. Recently Sollich [59] has determined several FDRs of quadratic 
operators (in the spin s{) that could be compared with FT results. He considered the 
bond energy observable B\ = |(si — Sj) 2 , the product observable Pi = SiSj (with i,j 
nearest neighbours) both in the real and in the momentum space, and the total energy 
E. The exact results of his analysis are that Xp° = X|? = X^ ^ X|? . The scaling 
forms for these observables have been also derived and they agree with FT expectations 
equations (106) and (107) with ap = — 1 — d/2, ap = 1 — d/2 and ap = d/2 — 3 in 
d < A (and ap = ap = as + 2 = 1 — d/2 in d > A). The calculation of response and 
correlation functions of high-order observables is rather cumbersome and it has not yet 
been carried out completely (see for details [23]). 

Two-dimensional Ising model: The physically relevant cases with finite N and 
d = 2, 3 are not analytically solvable and have been intensively investigated by means 
of Monte Carlo (MC) simulations. For the two-dimensional Ising model the dynamical 
critical exponents are known with very high precision from lattice techniques (we recall 
that the static model is solvable, see, e.g., [102]). The most accurate determination 
of the critical exponent z is z = 2.1667(5) [103] (this number is obtained by means of 
equilibrium MC, since the exponent z is the same in- and out-of-equilibrium) . This value 
is in agreement with other less precise simulations [104, 105, 106, 107], but it is somehow 
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bigger than the most precise FT estimates, i.e., the three-loop e expansion [108], four- 
loop fixed dimension [109, 110] and the interpolation of two-loop result between one and 
four dimensions [111]. The most accurate value of the non-equilibrium critical exponent 
9' has been obtained in reference [105], i.e., 9' = 0.191(3) that, using the scaling relation 
(95) leads to 9 = 0.383(3), compatible with other MC results [46, 106, 107] and again 
slightly bigger than two-loop FT estimate [24] . 

Concerning the FDR, earlier investigations indicate X°° = 0.26(1) [17]. Recent 
simulations provide results of increasing accuracy thanks to the improved methods that 
have been introduced to measures the quantities of interest. From the analysis of zero- 
momentum observables the very accurate estimate X°° = 0.340(5) has been obtained 
[58] (the method employed in [58] has been questioned, see for details [112, 113]). 
Recently a new algorithm has been proposed by Chatelain [114] and subsequently revised 
in [115]. In standard simulations some integrated form of the response function is 
measured after the application of a magnetic field. The new algorithm instead allows 
one to measure the integrated linear response function without applying any magnetic 
field to the system, avoiding in such a way an eventual crossover towards the non- 
linear regime. By applying this algorithm, the result X°° = 0.33(2) has been found 
[114]. The study of a more general model that belongs (apart one exceptional point 
in its phase diagram) to the Glauber-Ising universality class, even if it does not satisfy 
detailed balance, gives X°° = 0.33(1) [64]. A few attempts have been made to measure 
the critical FDR for the energy X^. In [64] the standard algorithm produced a signal 
too noisy to have any reliable result. By considering, instead, coherent observables 
X^ = 0.33(2) has been obtained [58]. This result led the authors of [58] to put forward 
the conjecture that X°° is the same for all the observables, according to the idea of 
the existence of a unique effective temperature, as it is the case for mean-field glassy 
models [5]. We comment that along the same lines of references [114, 115], an alternative 
algorithm measuring directly the response function has been proposed [116], but it has 
only been tested in one-dimensional systems. 

Three-dimensional Ising universality class: In this case the numerical simulations 
are obviously less accurate than in two dimensions. The values of the static critical 
exponents are reviewed in reference [67] and will be not reported here, neither will for 
the O(N) models with N > 1. The best determinations of the equilibrium dynamic 
critical exponent are z = 2.032(4) [105] and z = 2.055(10) [117]. They are not in 
agreement within three standard deviations, but they are anyway almost compatible 
with the result of less accurate simulations [118, 107]. Again the numerical values are 
a bit higher than the best FT estimates giving z ~ 2.02 [108, 109, 110, 111]. For the 
initial slip exponent the best value is 9' = 0.104(3) [105] in agreement with less precise 
estimates [46, 118] and with the two-loop e expansion [24, 25]. Through the scaling 
relation (95) and using the safe value z = 2.04(2), this value of 9' leads to 9 = 0.14(1). 
These values confirm that the e expansion converges better in d = 3 than in d = 2. 

Regarding the FDR, only the preliminary analysis of reference [17] is currently 
available, giving X°° ~ 0.40. 
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Two-dimensional XY universality class: The non-equilibrium dynamics of this 
universality class has been studied in [16, 119, 120, 121, 122, 123]. The model is quite 
peculiar since, despite the absence of a spontaneous magnetization, it is critical for all 
temperatures below the so-called Kosterlitz-Thouless temperature Tkt (i-e., in the RG 
language, it has a line of fixed points) with a continuously varying critical exponent 
r](T) that ranges from to 1/4 as T changes from to T KT [124]. Regarding the 
exponent z, all the approaches give z — 2, but with logarithmic corrections to the 
scaling depending upon the presence of vortices in the initial state [121]. For a quench 
from a high-temperature state to one of the critical points, numerical simulations have 
been so far inconclusive in determining non-equilibrium exponents. Indeed in [120] 
the value of the critical exponent 6' was found to vary continuously along the line of 
critical points, whereas in [123] it turned out to be temperature independent. A precise 
determination of this exponents surely deserves further analytical and numerical study. 
Even the determination of X°° is problematic: In fact a standard FD plot provides the 
approximate relation X°° = 1/2(T/T KT ) for T < T KT [123], whereas a direct evaluation 
gives X°° ~ 0, independently of T [123]. Given this argument has not yet been settled 
we do not report these results in table 1. Let us mention that the FDR has been studied 
also for the quench between two (critical) points below T KT [122, 125, 123]. 

It is worth noting that Schehr and Le Doussal [126, 60] obtained a temperature- 
dependent FDR (similar to that just mentioned) X°° = (2(1 + e 7B r) + O^ 2 )^ 1 , where 
t = (T g — T)/T g (being T g the temperature of the glass transition in the model) in the 
glassy phase of the two-dimensional Cardy-Ostlund model [127], that has a line of fixed 
points as well. 

Three-dimensional O(N) models with N > 1: These models have attracted a 
modest interest, mainly because the most relevant experimental realizations of O(N) 
models belong to different dynamic universality classes as discussed in section 3.4. 

To our knowledge the three-dimensional XY model with purely dissipative dynamics 
has been numerically studied only in [128], using the approach discussed in [114, 115]. 
X/z ~ 1.34 was obtained, without providing independent determinations of the two 
exponents. A direct determination of the FDR X(t, s) in the long-time limit turns out to 
be unreliable, since the numerical data are still affected by large statistical fluctuations. 
On the other hand, a linear extrapolation to the small-s/t region provides the estimate 
X°° = 0.43(4). 

Concerning the 0(3) universality class, the determination of z has been attempted 
only in [129], suggesting values slightly smaller than 2, i.e., z = 1.96(6) [129]. On the 
other hand, FT studies unambiguously give z > 2 [108, 130, 131]. In fact, we point out 
that the early three-loop FT computation of z [132] is affected by a numerical error, 
later corrected in [108], that led to the erroneous conclusion z < 2 in three dimensions. 
The correct value is z = 2 + Ri], with a value of R that is independent of N up to 
three loops and reads R = 0.726[1 — 0.1885e + 0(e 2 )] [108]. A more accurate numerical 
determination of z surely deserves further analysis. 
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Check of the universality ofX°°: The issue of universality of X°° has been addressed 
in [133] by using the algorithm of [114, 115] for a variety of models, some of which 
are not related to the O(N) models considered here. The long time non-equilibrium 
dynamics of the Ising, 3-state clock [133] and 4-state Potts [133, 102] models in two 
dimensions has been studied on square [sq], triangular [tr] and honeycomb [he] lattices. 
The values found for X°° are compatible with each other and with those already 
reported in the literature. For the Ising model it has been found X°°[sq] = 0.330(5), 
X°°[tr] = 0.326(7) and X°°[hc] = 0.330(8) on square, triangular and honeycomb lattice, 
respectively. For the 3-state clock model, instead, X°°[sq] = 0.406(6), X°°[tr] = 0.403(5) 
and X°°[hc] = 0.401(6), whereas for the 4-state Potts model X°°[sq] = 0.475(10), 
X°°[tr] = 0.470(12) and X°°[hc] = 0.467(7). These findings support the universality 
of X°° . Moreover, Chatelain measured X°° for other models on the square lattice, 
which at least in equilibrium are expected to belong to the same universality class as 
the previously mentioned models. Namely, he investigated the Ashkin- Teller model 
[133, 102] at the point of its phase diagram where it belongs to the same universality 
class as the 3-state clock model, finding X°° = 0.405(10), in perfect agreement with the 
previous results. Then he considered the Multispin Ising model [133] and the Baxter- Wu 
model [133, 102] (defined on the triangular lattice and belonging to the same universality 
class as the 4-state Potts model) finding X°° = 0.467(10) and X°° = 0.550(15), 
respectively (note that, to our knowledge, the latter value is the first evidence for 
X°° > 1/2 in ferromagnetic models quenched at their critical points from an initially 
uncorrelated state). Whereas the former value is in perfect agreement with the result for 
the 4-state Potts model, the latter is definitely incompatible with it. This can be a signal 
of a non-universality of X°° (according to [133], even the dynamical exponent 9 seems to 
depend on the specific lattice model, within the same dynamic universality class), or of 
the existence of a different (non-equilibrium) universality class to which the the Baxter- 
Wu model belongs, or simply of large cross-over effects that make particularly difficult 
to reach the asymptotic scaling behaviour. It is still unclear whether these possibilities 
can explain the (somehow) surprising results of [133]. Further investigations and cross- 
checks are surely required. 

4-1. Field- Theoretical Approach 

The coarse-graining of the lattice models specified by the Hamiltonian (110) with a 
Glauber-like dynamics leads to the purely dissipative dynamics of a iV-component field 
tp(x.,t) (Model A of [26]), described by the stochastic Langevin equation (33), where 
T~i[ip] is the ip A Landau-Ginzburg Hamiltonian given in equation (22). The equilibrium 
correlation functions, generated by the Langevin equation (33) and averaged over the 
noise (, can be obtained, as summarized in section 3.3, by means of the field-theoretical 
action (35). The problem is completely specified once equation (68) is assumed to 
provide the statistical distribution of the initial condition (with Tq 1 = given we are 
interested only in the leading scaling behaviour, and a(x) = 0). The propagators (i.e., 
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Table 1. Estimates of X°° in purely dissipative models quenched from a completely 
disordered initial state, [tr] and [he] stand for simulations on the triangular and 
honeycomb lattices, respectively. When not indicated, the simulation is on the square 
lattice. Models belonging to the same universality class are marked with t and 



Model 


Reference 


Method 


X°° 


Random Walk 


[161 


Exact 


1/2 


Gaussian model 


[161 


Exact 


1/2 


d-dim Spherical Model 


[171 
l x 1 j 


Exact 


1-2/d 


1-dim Tsirif? Model 


[95 941 


Exa ct 


1/2 


2-dim Ising Model 


[171 
l x 1 j 


MC 


0.26(1) 




[581 


MC 


0.340(5) 




[1141 


MC 


33 ( 2) 




[641 


MC 


0.33(1) 




[1331 


MC 


0.330(5) 




[1331 


MC[tr] 


0.326(7) 




[1331 


MCfhcl 

111 * — llv 


330(8) 




[201 


FT 0(e 2 ) 


0.30(5) 


O-A im p\ t p pi nplc A/ToHpIT 


fl 331 


MC 


4-06(6) 




[1331 


MC[tr] 


0.403(5) 




[1331 


MC[hc] 


0.401(6) 


2-dim Ashkin-Teller Model 1 ' 


[1331 


MC 


0.405(10) 


2-dim 4-state Potts Model* 


[1331 


MC 


0.475(10) 




[133] 


MC[tr] 


0.470(12) 




[133] 


MC[hc] 


0.467(7) 


2-dim Baxter- Wu Model* 


[133] 


MC[tr] 


0.550(15) 


2-dim Multispin Ising Model* 


[133] 


MC 


0.467(10) 


3-dim Ising Model 


[17] 


MC 


~ 0.40 




[20] 


FT 0(e 2 ) 


0.429(6) 


3-dim 0(2) Model 


[20] 


FT 0(e 2 ) 


0.416(8) 




[128] 


MC 


0.43(4) 


3-dim 0(3) Model 


[20] 


FT 0(e 2 ) 


0.405(10) 


3-dim Random Ising Model 


[21] 


FT Q(y/e) 


~ 0.416 



the Gaussian two-point correlation and response functions) are obtained in a standard 
way, by considering the quadratic part of the action: 

B°(t,s) = n(^(q,s)^(-q,t)) G = nS ij 9(t-s)e- n ^ +r ^ t - s \ (112) 

e -Q(c l 2 +ro)\t-s\ _ e -Q(ci 2 +r Q )(t+s) ^ 

where the subscript "G" is used to remind that we are considering the Gaussian 
approximation. The response function (112) is the same as in equilibrium. Equation 



= (v>i(q,s)¥>j(-q,*))G = 



a 2 4- rn 
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Table 2. Estimates of Xq for observables different from the magnetization. The 
meaning of the symbols [u], [1] and [q] in the FT estimates is explained in section 4.4 



Model 


Reference 


Observable 


Method 




Gaussian model 


[23] 


All 


Exact 


1/2 


(i-dim Spherical Model 


[59] 


E 


Exact 


^ 1 - 2/d 




[59] 


B,P 


Exact 


1-2/d 


1-dim Ising Model 


[58] 


E 


Exact 







[58] 


Other 


Exact 


1/2 


2-dim Ising Model 


[58] 


E 


MC 


0.33(2) 


2- dim Ising Model 


[23] 


E 


FT 


0.33[u], 0.20[1] 


3-dim Ising Model 


[23] 


E 


FT 


0.39[u], 0.37[1] 


3-dim 0(2) Model 


[23] 


E 


FT 


0.37[u], 0.31[1], 0.22[q] 




[23] 


T 


FT 


-0.41 [u] 


3-dim 0(3) Model 


[23] 


E 


FT 


0.35[u], 0.30[1], 0.21[q] 




[23] 


T 


FT 


-0.41 [u] 



(113) reduces, instead, to the equilibrium form when q ^ and both times t and s go 
to infinity while t — s is kept fixed. This ensures that quasi-equilibrium (in the sense 
discussed in section 2.5) holds for the order parameter, at least within the Gaussian 
approximation. 

4-2. Gaussian Fluctuation- Dissipation Ratios 

In the Gaussian model (equation (22) with go = 0) the response and correlations 
functions are exactly known and therefore the FDR can be easily evaluated. From 
equations (112), (113) and the definition (18) one finds 

R 1 

^ q(t ' S) = ^ = 1 + e -2n(q 2 +ro)* ' (114) 

where we set s < t and removed the subscript "G" since in this section we are only 
concerned with the Gaussian approximation. If the theory is non-critical (r > 0), the 
limit of the FDR for s — > oo is 1 for all values of q, in agreement with the idea that in 
the high-temperature phase all the fluctuating modes have a finite equilibration time. 
Equilibrium is hence recovered and the FDT applies. In the critical theory r = 0, if 
q 7^ 0, the limit ratio is again equal to 1, whereas for q = one has Xq=o(t, s) = |, 
independently of s and t. This shows that the only mode that "does not relax" to the 
equilibrium is the zero mode in the critical limit. This picture has been confirmed by a 
one-loop computation [19]. 

In order to understand why measuring observables in momentum space is more 
effective than doing the same in real space [58], it is instructive to look at the FDR 
in real space. Computing the Fourier transform of equations (112) and (113) one finds 
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(117) 



(with Q = 1) 

R*(t, s)= J (dg)e iqx e-( q2+rt ^*- s ) = j^y^{t ~ s )- rf / 2 e -«(^ s )- x2 /[4(t- s )] ? ( 115 ) 
C x (f,a) = /(dg)e^ — . (116) 

J q 2 + 

The derivative of the C x (t, s) with respect to s is 
<9 s C x (t, s) = ^ d/2 
Therefore the FDR is 

(. \ d/2 2 

^) e- 2 ''» s+ W^I. (118) 

Again, for fixed x and ro 7^ 0, in the ageing regime (t, s — > 00, fixed s/£) this quantity 
reaches the asymptotic value Xjf = 1, and the system equilibrates. Instead for r = 
and t, s — > 00 with fixed s/t it reduces to 

X x 1 (t, S ) = l+ (L-4 , (119) 



j + s, 

explicitly showing the crossover from the equilibrium value 1 at s/t — 1 to the pure 
ageing behaviour X°° = 1/2 for t ^> s (apart from d = where the model is a random 
walk) . 

The above considerations illustrate that XJf displays quite a long transient, lasting 
for a time of the order of the age of the system, before reaching its asymptotic value. 
This transient is absent in Af q=0 , resulting in a more effective determination of X°°. As 
we shall see, the corrections to Xq(t, s)/X°° due to the interaction term g are very small 
in magnitude (at least for Model A dynamics of O(N) models), making this remark still 
valid. 

For the Gaussian theory the FDR can easily be computed for a generic observable. 
In particular it was shown [23] that the FDR for a set of local one-point observables 
is always equal to 1/2 in the long-time limit. This argument is very simple and we 
briefly recall it here: Let us consider the operators of the form O i>n = d l (p n . The critical 
response and correlation functions of the order parameter in real space are given by 
equations (115) and (116). The two-point correlation function of Oq^ is given by the 
n — 1-loop diagram with the two points connected by n correlation lines (see [23] for 
details). In real space its expression is simply given by the product of the n correlators: 

C%(t,s) = c n [C x (t,s)} n , (120) 

(c n is the combinatorical factor associated with the diagram) whose derivative is 

d,C%(t, s) = c„n[C x (t, s^dsC^t, s) . (121) 

Analogously, the response function is obtained from that one contributing to the 
correlation function by replacing an order-parameter correlator with a response function: 

i£(t, s) = c n n[C x (t, sT-'R^t, s) , (122) 
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where the factor n comes from the response operator O , n = ^Co.n-i^ ( see equation 
(99)). Note that the combinatorial factor c n is the same as for the correlation function. 
The associated FDR is 

^•" s ja=TO sA - (M) - <i23) 

Therefore the FDR of powers of the order parameter in real space is equal, for all times, 
to the FDR of the order parameter itself. 

Before considering the effects of the derivatives let us consider the previous relation 
in momentum space. Remembering that the product of two functions after a Fourier 
transformation becomes a convolution, one has 

B%(t,s) = c n n(C*---*C*R) fl , (124) 
d s C°(t, s) = c n n(C * • • • * C * C") q , (125) 

where * is the convolution, C and R are in momentum space (with the time dependence 
understood), • • • means n — 1 times and C'^ = <9 s C q . For q = the previous relations 
become 

R% =0 (t, s) = c n n J(dp)(C*---* C) P R_ P , (126) 
d s C° =0 (t, s) = c n n J(dp)(C * • • • * C) p C'_ p . (127) 

Thus 

[y o M-i _ K d p)(c * ■ ■ ■ * c) p c_ p _ /(d P )(c * ■ ■ ■ * cyg_p^ 

Lq=ol ' JJ "/(dp)(C7*...*C7) p i2_ p - f(dp)(C*---*C) p itL p ' ^ 
i.e., the inverse of is a weighted average of XZ p , with weight (C * ■ ■ ■ * C) p i?_ p . 
In the limit t — > oo, s — > oo this weight is expected to be peaked around p = 0, giving 

= lim^oo lim^oo X® =0 (t, s) = lim^oo lim^oo X p =o(t, s) = X°°. Thus the relation 
between the different FDRs holds also when considering Co,n>i- Note that in momentum 
space, at variance with the relation (123) holding between FDR in real space, only the 
long-time limit of the FDR reproduces the FDR for the fields. 

Let us now account for the effect of the derivative. In momentum space these are 
simply multiplication by q* for each insertion, modifying the weight in equation (119) 
by a factor q 2 ' 1 which does not change our conclusion about the long-time limit. Since 
all correlation and response functions of all local operators can be written in terms of 
those of O it n, this concludes the argument. 

It is worth mentioning that for observables different from the order parameter, 
Afq (t, s) is expected to depend non-trivially on the times even in the Gaussian model. 
Hence the previous argument about the advantage of using observables in momentum 
space to determine X°° does not apply straightforwardly. 

4-3. Two-loop order parameter FDR 



The program of calculating the universal two-point functions of the order parameter, 
the associated non-universal amplitudes and the corresponding FDR has been worked 
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out up to two loops in an e expansion [19, 20]. The analytic manipulations required 
for the calculation are quite long and involved. We summarize here only the results, 
referring the reader interested in technical details to the original works [19, 20]. The 
obtained exponents z and 9 agree with the already known two- loop expressions [30, 24]. 
The scaling of the response function (see equation (86)) is characterized by the non- 
universal amplitude 



A R = 1 + e- 



3(iV + 2) 
8(iV + 8) : 



/(0)-4 7E log- 



+ 0(e 



and the universal function 



F R (v) = 1 + e 



2 3(iV + 2) 



[f(v)-f(0)] + O(e" 



(129) 



(130) 



8(iV + 8) 2 ' 

whereas the scaling function of the derivative with respect to s of the correlation function 
is (c.f. equation (88)) 

2 3(iV + 2) 



Fdc(v) 



~, 4, 1 - v T . . 

2 log -log h V(v) 

& 3 & 1+v K ' 



0{e*). 



(131) 



8(iV + 8) 2 

The correlation function has been computed only up to one loop [19]. We do not report 
here the quite long expression of Age that can be worked out from the results in [20] . In 
equations (129), (130) and (131) f(v) and ty(v) are regular functions for small v, defined 
for < v < 1. Explicit expressions can be found in [20]. For the present purposes it is 
enough to note that they are of order 1 up to quite large v, i.e., v ~ 0.9. In terms of 
Ar, Aqc and 9, the long-time limit of the critical FDR is obtained via equation (90): 

(X 00 )- 1 _ N + 2 9 N + 2 



1 + 



4(iV + 8 



-e + e z 



(N 



N + 2 3(3iV + 14) c 



8 



4(iV + 8) 



+ 0(e 3 ), (132) 



2 

with c = —0.0415 . . . (its analytic expression is given in [20]). 

Let us comment on these results and compare them with those currently available. 
The calculation shows that fluctuations give corrections to the Gaussian value of F R 
and Fqc of order e 2 . Unfortunately these are extremely tiny, making the numerical and 
experimental detection very hard. This is due to the factor e 2 [3(iV + 2)}/[8(N + 8) 2 ] in 
front of f(v) and W(v). For the Ising model, this factor is ~ 0.014 for d = 3 and ~ 0.05 
for d = 2, i.e., the two-loop corrections to the universal scaling functions are about 1% 
and 5% of the leading term in d = 3 and d = 2 respectively. This is not surprising given 
that even for the equilibrium scaling functions, the corrections to the Gaussian results 
involve the same prefactor [134]. 

However the tiny correction in Fr(v) is very important since it disagrees with the 
prediction, based on the theory Local Scale Invariance (LSI) [135], that F R (v) = 1 
exactly. The LSI is a proposed extension of conformal invariance to system with 
anisotropic scaling, e.g., to dynamics where time and space scales differently. Given 
the large success of conformal invariance (see, e.g., [57]) for isotropic scaling it is very 
interesting to provide generalizations to the anisotropic case. For Glauber-Ising model 
in d — 2, 3, the agreement of LSI predictions for the integrated response function with 
numerical simulation is really remarkable [54]. However, measuring a more suitable 
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integrated response function, a small deviation from LSI has been recently detected 
[136]. This deviation is in qualitative agreement with the FT prediction in equation 
(130). The failure of LSI has probably to be attributed to the limits of applicability of 
the theory [137, 125] and might be traced back to the fact that the term Q (p 1 (generated 
by a non- vanishing thermal noise) breaks the Galilei invariance of the action (35) [125]. 
However, it has been argued that LSI should give exact predictions whenever z = 2 
[125], even in the presence of the thermal noise. Field-theoretical results agree with this 
observation (see also sections 7 and 8). In fact, in the e-expansion, Fr(v) is Gaussian at 
least up to the order at which z = 2 (i.e., 0(e) in Model A of (f 4 theory, O(e ) in Model 
C, and 0(^/e > ) in Model A of the dilute Ising model). 

Let us remark that the time dependence of Xq=o(t, s) shows up only through the 
two-loop corrections appearing in equations (130) and (131). According to previous 
discussion, these are expected to be quantitatively small, making it difficult to measure 
the difference between Xq=o(t, s) and its asymptotic value X°°, as noted in [113]. 

The expression (132) of X°° for general N provides quantitative predictions for a 
large class of systems. To get some numerical estimates out of this two-loop expansion 
a direct summation (Pade approximant [2,0]) and an "inverse" one (Pade approximant 
[0, 2]) has been performed in [20]. From them some general trends may be understood: 

• decreasing the dimensionality, X°° always decreases (at least up to e = 2); 

• increasing N, X°° decreases, approaching the exact result of the spherical model; 

• for N = oo the curve of the [0, 2] approximant reproduces better the exact result 
in any dimension with respect to the [2, 0] approximant. 

The last point suggests the use of the [0,2] value as estimate of X°°, also for physical 
N. As an indicative error bar the difference between the two approximants has been 
used [20]. Accordingly, the numerical predictions for various cases can be summarized 
as follows: 

a) Ising model in two and three dimensions: Using the procedure just outlined, 
the FT estimates are X°° = 0.30(5) for the two-dimensional Ising model and 
X°° = 0.429(6) for the three-dimensional one [20]. In two dimensions the agreement 
is quite satisfactory despite the fact that a priori such a low-order perturbative 
expansion is not expected to provide such nice results. The estimate X°° = 0.30(5) 
is marginally compatible with the first numerical determination X°° = 0.26(1) 
[17]. The agreement with the subsequent estimates reviewed in table 1, and giving 
X°° ~ 0.33, is clearly improved. In three dimensions the preliminary investigation 
of [17] gives X°° ~ 0.40 for the Ising model. Albeit quite inaccurate, it is in nice 
agreement with the FT estimate X°° = 0.429(6). Further simulation of this model 
in three dimensions would provide an important test of the actual reliability of this 
quite precise FT estimate. 

b) XY model in three dimensions: The result for X°° yields also the accurate prediction 
X°° = 0.416(8) for the purely dissipative three-dimensional XY universality class. 
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Currently, its ageing properties have been investigated only in [128] where the value 
X°° = 0.43(4) is reported. On the one hand this result is compatible, within the 
estimated error bars, with the FT estimate. On the other hand, its numerical 
accuracy is still not sufficient to distinguish the differences with the Ising and 
Heisenberg universality classes in three dimensions. 

c) Heisenberg model in three dimensions: For the Model A Heisenberg universality 
class in three dimensions, the estimate of X°° provided in [20] is X°° = 0.405(10). 
Currently, no numerical results are available for this case. 

4-4- Two-loop FDRs of quadratic operators 

The two-loop response and correlation function of the energy and the tensor operator 
(see equations (103) and (104), respectively) were computed for generic times in reference 
[23] . Their general expression is very cumbersome and we refer the interested reader to 
the original reference. Here we report and discuss only the value of Xq [23]: 



Comparing these results with equation (132) one concludes that the long-time limit of 
the fluctuation-dissipation ratio depends on the particular observable chosen to compute 
it. 

Let us discuss in more detail some specific cases: 
• N = oo: In this case one finds 



which is clearly different from the expression for the FDR of the order parameter, 
as can be checked expanding — 1 — 2/d close to d — 4. We note that equation 
(134) agrees with the expansion close to d = 4 of the result for the total energy 
in the spherical model [59]. On the basis of combinatorial arguments one expects 
X%? to coincide with X^ for N = oo to all orders in e [23]. This expectation is 
confirmed by the two- loop FT computation equation (133). 

These findings agree with results obtained for the spherical model [59]. 
Unexpectedly, the observable P of the spherical model (see previous section) can 
not be naively identified with (p 2 for N = oo. On the other hand 
agrees with the naive identification of B as the Laplacian of P in the continuum 
limit. This fact calls for a more rigorous investigation of a possible correspondence 
between the two models, beyond the case of equilibrium dynamics. 

• Ising Model: The direct estimate (unconstrained) from the O(e) series in equation 
(133) for e = 2 (giving X|?[u] ~ 0.28) is probably unreliable, as was the case for 




(133) 




(134) 
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X™, i.e., X~[lloop] ~ 0.42 [19] whereas X^[21oops] = 0.30(5) [20]. To obtain a 
more reliable result without computing the 0(e 2 ) term, it was proposed in [23] to 
constrain linearly the available one-loop series to assume the exactly known value 
for d — 1 (i.e., e = 3). Assuming a smooth behaviour in e up to e = 3, this gives 
[23] 

A^[l] = i(l-£) [l + ^ + 0(e 2 )] , (135) 

which has the same e-expansion as equation (133), but it is expected to converge 
more rapidly to the correct result. From equation (135) for the two-dimensional 
Ising model one gets -X"|?[l] ~ 0.20, that is much lower than the value determined 
so far, — 0.33 [20, 114, 58, 64, 133]. However, a robust field-theoretical 
prediction of X|? for the two-dimensional Ising model requires a difficult higher- 
loop computation. For the three-dimensional Ising model equation (135) gives 
X|?[l] ~ 0.37, to be compared with the direct estimate X|?[u] ~ 0.39. Note that, 
as usual in d = 3, the spreading of the different estimates is much smaller, signaling 
a higher reliability of these predictions. 

O(N) model with N > 1: From equation (133) predictions can be obtained for 
the purely dissipative three-dimensional O(N) models with arbitrary N, again 
considering constrained analyses at the lower-critical dimension, which is d\ c d = 2 
in this case. For N > 2, X^(d = 2) = 0, since for d = 2 these systems are in the 
coarsening regime. In [23] this was assumed generically for N > 2. Within this 
assumption a linear constraint gives [23] 



A1?[l] = i(l-0 



16 -N 2 , 

1 + WTsf + 0{e) 



(136) 



On the other hand the exact results for the spherical model [59] suggest that the 
approach to d = 2 is quadratic rather than linear. For this reason, also the quadratic 
constraint [23] 

1 A e\2 



*f[q] , , ■ 2 



20 + AT 2n 

1+ 3wir +0(e) 



(137) 



was implemented. Comparing with the result for the spherical model, it was shown 
that for N = oo equation (137) provides a better approximation of the exact result 
than equation (136) . For the three-dimensional XY model (N = 2) one gets 
-X"J?[u] ~ 0.37 from direct estimate, -XJ?[1] ~ 0.31 from linear constraint, and 
X|?[q] ~ 0.22 from the quadratic one. For the three-dimensional Heisenberg model 
(N = 3) one finds, instead, Xf[u] ~ 0.35, Xf[\] ~ 0.30 and Xf[q] ~ 0.21. 
Even the results with constraints are rather scattered, making it difficult to provide 
firm estimates in d = 3. However, the conclusion X|? < X^ for all 2 < d < 4 is 
quite robust. Furthermore the difference between the two FDRs should be large 
enough to be observed in three-dimensional Monte Carlo simulations. The analysis 
of the non-equilibrium behaviour within the e = d — 2 expansion [30, 131] may 
clarify which, between the linear and the quadratic constraint, is the proper one 
close to d — 2. 
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We note that X™ is very close to making the numerical detection of such a 
difference probably very difficult. 

Let us contrast the FT results of [23] with those obtained for the one-dimensional 
Ising model [58, 97] and for the spherical model [59], where Xq = X^ for all the 
observables O, except for the total energy On the other hand, one could doubt that 
the energy operator is not as suited as others to define the effective temperature, being 
conjugated to the temperature of the bath but not to the actual one (if any) of the 
system. Nevertheless FT results show that there is at least one further operator, namely 
Ty, having X£P ^ X^, X^?. This supports the idea that a unique effective temperature 
can not be defined for this class of models [23] . 

We recall that in disagreement with FT calculation, the MC simulation of the two- 
dimensional Ising model apparently gives X^ = X^ [58]. Probably, a more accurate 
measure of X|? is required to detect the difference (if any in d — 2) between X%? and 

5. Surface behaviour 

It is a well-known experimental fact and theoretical result that the presence of surfaces 
influences both the static and dynamic equilibrium behaviour of critical systems 
[80, 138, 139, 140, 141, 142, 143]. Upon approaching the critical point, local quantities 
such as the order parameter may display algebraic singularities at the surface differing 
from those characteristic of the bulk behaviour. Accordingly, one usually introduces 
surface critical exponents [80] in addition to the standard bulk ones. In the generic 
case, the former can not be expressed in terms of the latter. 

5.1. The semi-infinite O(N) model: statics and dynamics 

As an illustrative example, let us consider the influence of a planar surface on the model 
(110) defined on a semi-infinite hypercubic lattice. Its Hamiltonian is given by 



where r > is the coupling between nearest-neighbouring spins located at the surface. 
Because of the missing bonds, these spins experience an effective local magnetic field 
which is smaller than in the bulk, at least for sufficiently small r. Accordingly, 
one expects that the surface is not able to order independently of the bulk. Upon 
approaching the bulk critical point at the critical temperature T c , the surface and the 
bulk order simultaneously. This defines the so-called ordinary transition (OT) and the 
associated surface critical behaviour is characterized by a given set of surface critical 
exponents [80]. On the other hand, upon increasing r at fixed T > T c , it is possible to 
overcompensate the missing bonds at the surface. As a consequence, the surface may 
order in the presence of a disordered bulk, provided that d — 1 is larger than the lower 
critical dimension of the model. For r greater than the threshold value r sp there is a 




(138) 



(ij)ebuik 
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finite surface transition temperature T s (r) > T c at which the so-called surface transition 
takes place. The surface critical behaviour of the surface transition is the same as the 
bulk critical one of the d — 1-dimensional system. For T s (r) > T > T c the system is 
characterized by an ordered surface and a disordered bulk. Upon further decreasing 
T, the bulk orders at T c in the presence of an already ordered surface. This marks the 
so-called extraordinary transition. In the (T, r)-plane the three lines of phase transitions 
meet at the so-called multicritical special transition (SpT). 

The same scenario is valid in equilibrium and non-equilibrium dynamics in semi- 
infinite systems: The dynamical behaviour is affected by the presence of the surface and 
depends on the surface universality class. On the other hand no new and independent 
exponents show up in addition to the static surface and static and dynamic bulk ones. 
This is true in the case of both equilibrium [141, 142, 143] and non-equilibrium evolution 
[144, 145, 146]. 

The influence of a surface on the dynamical behaviour following a quench from the 
high-temperature phase to the critical point has been investigated in [147]. In particular 
the surface autocorrelation function and the associated surface autocorrelation exponent 
A C; i are determined analytically at the OT and SpT for the Model A dynamics of the 
O(N) Model up to the first order in the e-expansion and exactly for the O(N) model 
with N — > oo. MC simulations are instead used to investigate the same quantities 
at the OT of the two-dimensional Ising model with Glauber dynamics. This problem 
was reconsidered in [145], showing that the surface autocorrelation exponent is not 
a new dynamical quantity as claimed in [147] but it can be expressed in terms of 
the bulk autocorrelation exponent A c and surface and bulk static critical exponents. 
The two-time surface autocorrelation function C±(t,s) was studied in [146]. By using 
scaling arguments, it was shown that a novel short-time behaviour (referred to as 
"cluster dissolution") is expected in some cases, among which the three-dimensional 
Ising model at OT. In this case a non-algebraic decay of correlations is expected: 
\nCi(t, s) ~ —(t — s) K for t — s <C s where the exponent k is given in terms of already 
known static surface and dynamic bulk exponents. This scenario was confirmed by MC 
simulation [146]. 

Previous studies have been extended to the response function and to the full 
dynamic scaling regime in [148]. The effects of the surface on the ageing properties of a 
critical systems are investigated via MC simulations at the OT of the two-dimensional 
Ising model and at the OT and SpT of the three-dimensional one. In addition the 
Hilhorst-van Leeuwen model is also considered [148]. The scaling behaviour observed 
in all these models agrees with scaling arguments, both for the surface autocorrelation 
function Ci(t,s) and for the integral of the surface autoresponse function Ri(t,s). In 
terms of R\ and C\ the surface FDR can be defined as: 

TR x {t,s) 

Xl(M) = ^M)- (139) 
Xi is analogous to the FDRs previously introduced. As in the bulk, it is possible to show 
that Xi(t, s) is a universal function and its long-time limit X^° is a universal amplitude 
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ratio. As such, it is expected to depend on the specific surface universality class. In 
[148] it was found that = 0.37(1) for the two-dimensional Ising model (OT) whereas 
for the three-dimensional Ising model = 0.59(2) at OT and = 0.44(2) at SpT. 

From these results it is quite clear that X^° ^ X°° (c.f. table 1) and in particular 
that Xi depends on the surface universality class, as expected. 

It is worth noting that MC simulations suggest that X™ can be greater than \ at 
the OT point, as a difference with the bulk, where X°° is always < \ for O(N) models. 

5.2. The semi-infinite O(N) Universality Class: Field- Theoretical Approach 

In the previous sections we have discussed in detail how to apply the field-theoretical 
approach to study ageing phenomena in the bulk of a critical system. Taking advantage 
of previous works on the subject [147, 145, 148] we outline here a scaling approach to 
the non-equilibrium dynamics in the presence of spatial surfaces. Let us go back to 
equations (84) and (85), valid for a critical system in the bulk, i.e., far enough from 
any surface. From these equations it is possible to compute the two-point correlation 
and response functions in real space, which depend on the coordinates of the two points 
(x',t') and (x, t). In the following we assume t > t'. Because of the translational 
invariance in the bulk, this dependence is only on x' — x. Let us introduce in the system 
a planar d — 1-dimensional surface S. As a consequence, the translational invariance 
is broken in the direction perpendicular to S. For later convenience we decompose the 
position vector as x = (x||,xj_) where xy is the d — 1-dimensional component parallel 
to S and x± is the one-dimensional one perpendicular to it. The residual translational 
symmetry along S implies that the two-point correlation functions depend on Xy and 
xy through xj, — xy whereas they depend on x' ± and x± separately. The dependence on 
x' ± — x± is recovered only in the bulk, i.e., far away from S. For qy = (where qy is 
the momentum conjugate to xj, — xy) one can write equations (84) and (85) as 



R^= (t,t';x ± ,x' ± ) = (t-tr- 1/z (t/t') 6 F^x ± (t-t')- 1/z ,x' ± (t-t')- 1/z ,t'/t) , (140) 
C Q|=0 (t, f ; x ± , x' ± ) = t'(t - tr- 1/z (t/t') e F^x ± (t - t'y 1/z , x' ± (t - tT 1/z , t'/t) , (141) 



where a = (2 — r\ — z)/z. The superscript b reminds that these forms are valid, a 
priori, in the bulk. As discussed in [145] no genuinely new divergences arise when 
discussing the non-equilibrium dynamics in the presence of a surface. This means 
that the short-time expansion t' — > of the fields y?(x', t') and y5(x', t') (see equations 
(77) and (78)) is not influenced by the presence of S and, conversely, that their 
short-distance expansion for x' ± — > is not influenced by the presence of the time 
surface. Accordingly, the short-distance expansions (SDE) of </?(x',t') and y5(x',t') are 
the same (at least in Model A) as in the static case and have a SDE coefficient C*(x' ± ) 
[</?(x',f) ~ C*(x' ± )ip(x.' = (xj,^ = 0),t') and analogous for (/5(x',t')] that in the scaling 
regime is given by [80] 



(142) 
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The exponent (3i is defined analogously to the bulk exponent (3 and describes the 
behaviour of the surface magnetization mi upon approaching the critical point. For 
the O(N) universality class field-theoretical estimates of f3\ have been provided using 
different approaches. For recent reviews see [138, 139]. Equation (142) implies that 
for x± — > 0, -R q||= o, C q||= o ~ [x±Y /3l ~ /3 ^ u i.e., they have no analytic behaviour in this 
limit. The previous scaling forms can be written in a way that singles out this limiting 
behaviour. As a result the scaling forms in the presence of a surface are given by 
(hereafter we omit the subscript qy = 0) 

R(t,t';x ± ,x' ± ) = A S R (t-s) ai (t/t') e [xlx'J^-W" 

x F° R (x ± (t-t , )- 1/ *,x , ± (t-1?)- 1/z ,t , /t) , (143) 
C(t,t';x ± ,x' ± ) = A s c t'(t - t') ai (t/t') e [ Xl _x'^- p)/v 

x F^x ± (t-t'r 1/z ,x' ± (t-t')- 1/z ,t'/t) , (144) 

where a± = a — 1/z — 2(/3i — /3)/(vz). The non-universal amplitudes A S RC are fixed 
requiring F RC (0,0,0) = 1. The previous equations provide an extension of the scaling 
forms presented in [145, 148]. As a consequence of equation (144) one gets 

d t ,C(t,t';x ± ,x' ± ) = Al c (t -tT(t/t') e [xLx'jCfc-fl/" 

x F s dC (x ± (t-tr 1/z ,x' ± (t-tr 1/z ,t'/t) , (145) 

where the non-universal amplitude A S 9C is fixed requiring Fg C (0, 0, 0) = 1. It is easy to 
show that A s dC = (1-9)A S C . 

It is natural to define the FDR in analogy to what has been done in the bulk: 

R(t, t'] x±,x' ± ) 
d s C(t,t' 

Using the previous scaling forms <^q| |= o can be written as 

' ' L> Al c Fl c (xJt-t')- 1 '',^Jt-t')-"',t'/t) ■ [M '> 
This explicitly shows that X^ =0 , just like Aq = o, is a universal function, being the ratio 
of two quantities that have the same scaling dimensions. The specific form of <^qn=o i s 
expected to depend on the universality class of the transition at the surface. Moreover, 
in the long-time limit t >> t' and sufficiently close to the surface x±,x' ± <C (t — t') 1 ^ 
the leading asymptotic value X s ' 00 is approached independently of the distance from the 
surface, i.e., 

x "° s ,£^ ± x «,=° (t - « = ^ = (T^r • (148) 

This is a consequence of the fact that at criticality the influence of the surface deeply 
penetrates into the bulk as time goes by, whereas a bulk-like behaviour can be observed 
only for x±,x' ± ^> (t — t') 1 ^. Following the argument outlined in section 2.4 one can 
argue that the asymptotic limit of the FDR just introduced is the same as that obtained 
by considering in equation (146) the surface two-point autoresponse and autocorrelation 
functions, i.e., quantities in real space. On the other hand, as its bulk counterpart Xq =0 , 



K=o(t, t'; x ± , x' ± ) ee ^Sl t' t\ ■ (146) 
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the expression in the mixed representation (q»,x±) is more suited to be computed 
within the field-theoretical approach and to detect in numerical simulations the ageing 
behaviour associated with the slow dynamics of the slow mode qy = 0. 

Let us check the scaling relations (143) and (144) within the Gaussian model and 
compute the associated FDR (147) in the case of OT and SpT. In the presence of the 
surface the equilibrium correlation and response function are given by [80] 

At) = R^'^ (qy , x± — x' ± , At) ± # (b ' e) (qy, x± + x' ± , At) (149) 
C< s ' e > (qy , x ± , x' ± , At) = (qy , x ± - x' ± , At) ± (q,| , x ± + x' ± , At) (150) 

for S Q P J. RM (q\\,x±,At) and C (b ' e) (q||, x±, At) are the bulk equilibrium response and 
correlation functions in the (qii , xjj-representation and At = t — t' . Let us recall that 
in momentum space (we set 0,-1) they read 

i? (b ' c) (q, At) = e(At) e -^ +ro)At , (151) 

C (b ' e) (q, At) = — ^ e -( q2+r °)l A *l (152) 

q 2 + r 

[q= (<l\\,Q±)}- Accordingly, 

R^ c \q h Ax ± ,At) = [ ^R^ e \q,At)e iq±Ax± 

J 27T 

= ^(At)(4 7 rAt)- 1 / 2 exp |-(qj + r )At - (153) 

POO 

C (b ' e) (q||,Ax ± ,At) = / duR^ c) (q h Ax ± ,u) . [FDT] (154) 

J\At\ 

As reviewed in section 3 the effect of the quench from the high-temperature disordered 
phase amounts to the introduction of a time surface in the field-theoretical description 
of the model. At the fixed point r = oo the boundary condition at the time surface 
is a Dirichlet one. Accordingly, the correlation function is modified, compared to the 
time-bulk one, by the introduction of an "image term" in time, i.e., 

(q;M ) — 

C( b ' e \q,t — t') — C( b ' e )(q, t + t'), whereas the response function is not modified: 
(q; t, t') = i?( b ' e )(q, t—t'). The presence of a spatial surface in the problem does not 
change the previous considerations. Accordingly, one finds that for the non-equilibrium 
evolution (ne) after a quench 

R {s ' ne) (q.\\,x ± ,x' ± ;t,t') = R {s ' c) (q h x ± ,x' ± ,t-t') , (155) 
C^ ne \q\\,x ± ,x' ± ;t,t') =C( s ' e \q\\,x ± ,x' ± ,t-t')-C( s ' e \q\\,x ± ,x' ± ,t + t') . (156) 

By means of equations (149), (150), (153) and (150) it is simple to find that _R( s ' ne ) 
and C^' 116 ** have the expected scaling forms (143) and (144) with the proper (mean- 
field) values of the critical exponents. Moreover, the non-universal amplitudes are 
easily computed: A s R [SpT] = tt" 1 / 2 , A s R [OT] = (An)- 1 / 2 , A s c [SpT] = 27T- 1 / 2 and 
A s c [OT] = it- 1 ' 2 . Using equation (148) one concludes that the surface FDR is equal to 
\ within the Gaussian theory for both SpT and OT, as it is the case for the Gaussian 
bulk FDR. 
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The previous equations can be rewritten in real space. In particular, the FDT 
together with equation (156) yields 

d t ,C (s ' ne) (x,x';M') = i2 (B,e) (x,x',i-0+i2 (s,e) ( x . x '.* + • ( 157 ) 
Accordingly, the FDR in real space can be written as (t > t') 

[X s (x, x'; t, t')]- 1 = 1 + „. x ' x ' t + t (158) 

where i?( s - e )(x, x', At) is given by the Fourier transform of equation (149). Let us 
investigate the cross-over between the bulk-like FDR and the surface-like one. For 
simplicity we consider the case x = x'. Then it is easy to find that 

[X(x,x,M)j -l+^—J l±exp{-xi/[t-f]}- (159) 

For ^> t + 1' one recovers the bulk value of the FDR equation (128) both for SpT 
and OT. In the opposite limit x 2 ± <^t — t' for the OT one finds 

/f_ft \ d/2+l 

[X>(x,x ] tJ)]- 1 = l+{—) (160) 

that differs from the corresponding bulk quantity and displays the same cross-over from 
the quasi-equilibrium regime t > t' to the ageing regime t ^> t' . On the other hand, in 
the case of SpT, one finds the same expression as in the bulk, given in equation (128). 
These results for the surface FDR provide the Gaussian theoretical predictions for the 
quantity X 1 numerically investigated in [148]. In particular, as far as the asymptotic 
value is concerned, the Gaussian result is = |, both for OT and SpT. At first 
sight, the non-monotonic behaviour of X^°[OT] as a function of the space dimensionality, 
found in [148] , can be hardly accounted for within the field-theoretical approach. Further 
analytical and numerical investigations are required to clarify this issue. 



6. Dynamics of a conserved order parameter (Model B) 

In this section we discuss the effect of the conservation law of the order parameter. As 
explained in section 3.4 this defines the Model B universality class. Compared to the 
Model A, Model B is much less studied in the literature, mainly because, due to the 
conservation law, its critical slowing-down is much more severe and it is quite difficult 
to reach the asymptotic behaviour in numerical and experimental investigations. 



6.1. Review of numerical and analytic results 

For Model B dynamics, the critical exponent z is z = 4 — rj exactly [26] and it has 
been argued that A = d [149], as found in several numerical studies for the one- and 
two-dimensional Ising model [149, 150, 151, 152, 153]. 

Earlier claims that the one- dimensional FD plot (and hence the FDR) of Kawasaki- 
Ising model reduces asymptotically to the Glauber one [52], have been recently revisited 
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in [152] and the apparent equivalence between the two models has been traced back to 
the time range on which the analysis of [52] was based. 

The scaling form of the autocorrelation function has attracted particular interest. 
For t, s 1 it has been found to scale like [152] 

C x=0 (t,s)~r d / z C(ts- 4 '), (161) 

where C(x) approaches a non-zero constant for x — > oo whereas C(x) ~ x (d~^')/z f or 
i « 1 and = 1 + (2 — 77) / (A' — d). This scaling implies the existence of a new 
relevant length scale ~ s^l z which controls to the crossover between the two regimes. 
Based on MC simulations of the one- and two-dimensional Ising model with Kawasaki 
dynamics it has been conjectured that A' = d + 3/2 [152]. This value has been later 
revised by Sire, who calculated exactly the correlation function for the O(oo) model, 
finding A' = d + 2 [153]. By means of simulations of the Ising model in d = 1,2 with an 
accelerated dynamics (in the same dynamic universality class), the relation A' = d + 2 
has been checked and it has been conjectured to generally hold for all systems with 
conserved order parameter. 

So far numerical simulations, both for d = 1 and 2, have not been able to 
explore the long-time regime from which X°° is determined. However the conclusion 
^Kawasaki > -^Glauber f° r the critical FDR of the magnetization seems quite robust, at 
least in two dimensions [152]. This fact is somehow surprising if one interprets T/X°° 
as an effective temperature. Indeed one would expect that a slower dynamics should 
result in an higher T cff . 

6.2. Field-theoretical approach 

The Field-theoretical approach to Model B dynamics starts from the effective action 
(39). The presence of the gradient term (that in momentum space gives a factor q 2 ) 
makes the difference with Model A dynamics: In fact it leads to an interaction vertex of 
the form q 2 g that reduces the primitive divergences of all Feynman diagrams so that no 
new divergences are generated by the dynamics [30, 24]. Hence the critical exponents 
are exactly obtained on the sole basis of power-counting: fj = —rj (i.e., z = 4 — rj 
[30, 28]) and i]q = 9 = [24]. Note that the response function in this case is not given 
by (y5_ q (s)<^ q (t)) = G q (t, s), but rather by i? q (t, s) = aq 2 G^(t, s) (see equation (32)). 

6. 3. Scaling forms 

The general scaling forms may be obtained solving RG equations with the methods of 
characteristic functions, as outlined in section 3.5.2 for Model A dynamics. In particular 
equations (72) and (73) for the correlation functions are unchanged and they provide 
quite stringent predictions since one exactly knows that r) = —rj and r/o = 0. 

The scaling of the autoresponse function at criticality is simply get from these 
equations with 6(1, 1, 0) = d: 

R*=o(t, s) = -a V^v^x, t)<p(0, s)) | x=0 = A R (t - s)-^ z ^ R (s/t) ,(162) 
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where, as usual, one singles out the non-universal constant Ar, so that Fr{v) is a 
universal function with .Fr(O) = 1. 

The autocorrelation function at the fixed point Tq 1 = is obtained from 5(2, 0, 0) = 
d-2 + rj = d + 2 + z: 

C x=0 (t, s) = A c s(t - s)-^/ z ^ c (s/t) , (163) 

where the proportionality to s is explicitly shown to stress the fact that the Dirichlet 
boundary condition in time enforces C x (t, 0) = for Tq 1 = 0. As emphasized in section 
3.6, the long-time limit of the correlation function for finite s can be influenced by the 
non-vanishing initial correlations, i.e., by a finite To- From equation (94) we get, in the 
long-time limit and neglecting further corrections coming from r^ 1 

Ait) = C x=Q (t,0) = ro^MM^Cx)) c~ rrH- 5 ^* = T~H~ d/z , (164) 

where we used 5(1,0,1) = d from equation (73). Accordingly, for t ^> s, C x= o(t, 0) — 
r ~ 1 t~ d / z , whereas, from equation (163), C x=0 (t, s) ~ st~( d+2 ^ z . As a fundamental 
difference with Model A dynamics, these two terms do not scale in the same way with 
t for t ^> s and in particular the former is the leading one for large t. Accordingly 
the autocorrelation exponent (see equation (96)) is given by A = d, in agreement with 
what has been argued in [149]. Furthermore a cross-over behaviour is expected for 
finite s when the contribution coming from finite initial correlations is of the same 
order of magnitude as the contribution with Tq 1 = 0. Comparing the previously 
given asymptotic behaviours one concludes that the former dominates when T~ x t~ d l z ^> 
s ^-(d+2)/z ^ j e ^ w h en To st~ 2 / z <^ 1. This agrees with the prediction (161), with = z/2 
(equivalently A' = d + 2), providing a FT proof of Sire's conjecture [153]. 

From the scaling of the autoresponse function (see the definition in equation (98)) 
we instead get Ar = d + 2 ^ A. However, such inequality does not prevent the definition 
of a universal FDR, as one might naively think on the basis of the fact that this leads 
to a different scaling dimensions for R and d s C. Indeed the contribution to d s C coming 
from finite To, leading in the large-t limit of C x= o — T ~ 1 t~ d / Z , is expected to scale as 
T 1 t~ d l z ~ 1 , turning out to be subleading in d s C x=0 (t, s) ~ t^( d + 2 )/ z . In fact the derivative 
is 

d s C x=0 (t, s) = A dC (t - sp d+ ^ z T dc (s/t) , (165) 
where, being 9 = 0, Aqc = Ac- It is possible to define the universal amplitude ratio 

x °° = iim iim =4^=^- ( i6e ) 

s^oot^oo o s U x=0 (t, s) Aqc -Ac 
6.4- Gaussian results 

The Gaussian response and correlation functions read for t > s 

R°(t,s) = aq 2 <^(q, S M-q,*)> G = aq 2 e~^ 2+r ^~ s \ (167) 
Cj(t,s) = (f(q,s)<f{-q,t)) G 

„-<Tq 2 (q 2 +r ( ,)(t-s) _ .- CT q 2 (q 2 +r )(t+s) „ n 
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where we keep the correction coming from non-zero r 1 . The Fourier transforms of the 
critical functions for x = are (setting a — 1) 

itWM) = r M (t - sr (d+2)/4 , (169) 

C x=0 (t, s) = r_ M a(f - s)- (d+2)/4 /G (|) + ro" 1 ^ + s)~ d/4 , (170) 

where r M = r((d + 2n)/4)/[2(47r d / 2 )r(d/2)] and 

/ oW = (i-,)<Wi-»>™ / '-fl+')'^ / ', (m) 

which is regular for small x: /g(0) = (2 — d)/2. These forms agree with those obtained 
from the solution of the RG equations (see equations (91) and (92)) with z — 4. 

From the scaling derived in the previous subsection, we can define the universal 
FDRs in real and reciprocal space as in equations (15) and (18). The argument presented 
in section 2.4 is easily extended to the present case, concluding X£L = X£L> . Here 
we use the notation given that y? q =o does not change in time, being a conserved 

quantity. Accordingly, i? q =o and d s Cq=o vanish whereas the limit q — > of their ratio 
is well-defined and finite. On the other hand, in a finite system of linear dimension L, 
the smallest allowed mode is q m i n with |q m i n | oc 1/L. The limit q — > is explored by 
studying the mode with q = q min in the limit of large L. In any case to stay within 
the non-equilibrium regime one has always to consider times much smaller than L z , 
otherwise the system equilibrates. 

The Gaussian FDR can be readily computed using the previous expressions, finding 

X~\t, S) = 1 + (l - S!±I^ e -2«r>(q*+ro)._ (172) 

As for Model A, all the modes with q ^ relax to equilibrium, i.e., X^l = 1- At 
variance with the case of Model A, Xq_+ (t, s) reaches a non trivial (but not universal) 
value 1/(2 + ror^ 1 ) also for r^ ^ (non-critical case). The conservation law allows 
arbitrarily slow modes which in the case of a purely dissipative dynamics is present only 
at criticality. The universal critical zero- momentum FDR is again |. 

Since it was instructive for the Model A dynamics, let us give a look at the Gaussian 
FDR in real space. From the Gaussian response and correlation functions reported above 
one gets 

,._ x (d+2)/4 , x (d+2)/4 

X;lo(t, s) = l+ (j-^ s ) - r^B d (L-^ J (t + s)-^ , (173) 

where B d = Y(d/A + l)/r(d/4 + 1/2). The term proportional to t 1 is only a correction 
to the asymptotic form for t > s — > oo. Anyway, as in the case of Model A, X~i has 
a slower approach to the asymptotic value X°° = 1/2 than that displayed by X^ . 
This faster approach suggests the use of coherent observables to determine X°° also in 
the case with conserved order parameter. Although this involves the limiting procedure 
described above, which is expected to be more noisy than working in real space, we 
think that the absence of a long transient allows for a more efficient determination of 
X°°. 
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6. 5. Exact solution for N = oo 

The next step in understanding general properties of X°° is the exact calculation for 
N = oo. Some quantities were already considered in the literature [84, 153]. In this limit 
the interaction can be self-consistently decoupled as g (p 4 (x,t) — > g /nC x= o(M)y? 2 (x, t) 
as for Model A dynamics. At criticality this leads to the the response function 
i? q (t, s) = q 2 G q (t, s) where 

G q (t, s) = e -q 4 (*--) e 2^q 2 (^-v^) _ ( 174 ) 

The constant is determined by a consistency condition (see [84], section IV). Hence, 
from the definition of the correlation function, one gets [153] 

C q (t, s) = 2q 2 f duG^t, u)G q (s, u) + T^G^f , 0)G q (s, 0) , (175) 
Jo 

explicitly showing that all the exponents are the mean-field ones, as expected. The FDR 
Xq-+o(t,s) turns out to be 1/2 for all times and without corrections, as an important 
difference with model A dynamics, where X°° = 1 — 2/d for 2 < d < 4 and N = oo. 

6.6. FDR beyond the Gaussian approximation 

It is easy to realize that the one-loop diagrams do not contribute to the critical FDR 
for q — > 0. In fact, as in Model A, they are given by tadpoles [19]. Their one-particle 
irreducible parts are independent of the external momentum q. On the other hand, the 
vertex carries a factor q 2 . Therefore, in the limit q 2 — > 0, the contributions of these 
diagrams vanish. Possibly non-vanishing corrections (if any) could come from diagrams 
whose one-particle irreducible part depends on q (as the two-loop sunset diagram). So 
one concludes that 

X°° = i + 0(e 2 ). (176) 
2 

A more detailed analysis of this problem is in progress [154]. 

Let us compare this result for Ag° wasaki with that one for Model A (AoJ auber ). By 
using equation (132) one finds that X- awasaki /X- aubcr = 1 + (N + 2)/[4(N + 8)}e + 0(e 2 ). 
Therefore, for small e > one has Ag° wasaki > A^ auber a relation that should not be 
changed by higher-order terms up to rather large e. This expectation is in agreement 
with the numerical results of [152], where such a relation has been observed for the 
two-dimensional Ising model. 

7. Coupling of a conserved density to the non-conserved order parameter 
(Model C) 

We now address the problem of when and how the critical behaviour is influenced by 
coupling a non conserved order parameter to a conserved density. As discussed in 
section 3.4 such a coupling characterizes the so-called Model C universality class. The 
starting point is the effective dynamic action (45) defined in terms of the two fields tp 
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and e. The propagators (Gaussian two point correlation and response functions) of the 
resulting theory are the same as Model A (equations (112) and (113)) as far as the order 
parameter (and the associated response field) are concerned, whereas, for e and e one 
has [81] 

(e(q, s)e(-q, t)) G = R° q (t, s) = 6{t - s)G £ {t - s), (177) 
(e(q, s)e(-q, t)) G = C° g (t, s) = G £ (\t - s\) + (c - l)G £ (t + s), (178) 

with G £ (t) = e- pQq2t . 

Let us briefly recall the scenario of fixed points for non-equilibrium Model 
C [155, 156, 81]. The fixed-point values for the couplings g and 7 are determined 
only by the statics. One has g* = g* A + Q^ 2 *, where g* A = 6e/(N + 8) + 0(e 2 ) is the 
fixed-point value of the coupling constant for Model A [30]. 

The value of 7 at the stable infrared fixed point depends on the sign of the specific- 
heat exponent a: 

_ 2it j , stable for a < , case I , 

[ Ar(iv+8) e + ^( e2 ) ' stable for a > , case II , 

in the case I, the dynamics of the conserved density decouples from that of the order 
parameter and one gets back to Model A (at least asymptotically). At the leading order 
in e-expansion one has, for the 0(N) model [30], 

a = WT8) t + 0i€2) - (180) 

therefore the truly Model C dynamics fixed point is stable for N < 4 + 0(e). In three 
dimensions, theoretical and experimental investigations showed that a is negative for 
N > 2 (see [67] for a comprehensive review). As a consequence the Model C dynamics 
may be realized only for the three-dimensional Ising model (N — 1) which has a positive 
a [67]. In two-dimensions, the boundary value between positive and negative a is TV = 1, 
so no physically relevant genuine Model C dynamics exists in this case. Albeit for d = 3 
with N > 2 and d = 2 with arbitrary N the dynamics is asymptotically Model A, a 
strong crossover is expected, especially in models with small a [156]. 

As far as p is concerned we have two possible stable fixed points determined by the 
equilibrium dynamics [156] 

(a) p* = 00, stable for N > Aq(e) = 4 - [15/4 + 3/2 log(4/3)]e + 0(e 2 ); 

(b) p* = 2/N - 1 + 0(e), stable for N < Aq(e). 

A third region with p* = 0, found in earlier calculation [155], has been recently proved 
to be only an artifact of the e-expansion and consequently it does not exist in any 
dimension [156]. This (fake) fixed point with p* = has been mentioned and discussed 
in many papers so far [26, 30, 81, 22]. Any reference to it can be ignored. 

Finally, as far as the non-equilibrium dynamics is concerned, it has been shown 
that, whenever a > 0, the fixed point value for c is c* = [81]. 
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We focus here our attention on the only relevant stable fixed point of the model, i.e., 
11(b). At this fixed point z = 2 + ajv exactly [30], so that z can be obtained from the 
precise estimates of static exponents [67] . 9 instead is an independent exponent known 
up to two loops [81] in e expansion. In [22] we obtained the one- loop non-equilibrium 
universal scaling function of the response 

Fr{v) = 1 + 4{N ^ N _ 1) Ml 1 + l)v) N - 2 (l — (N — l)v) N ]e + 0(e 2 ) , (181) 

that renders, for the physically relevant case of iV = 1, 

F R {v) = l-e V - + 0{t 2 ). (182) 
6 

Also the scaling function F c {v) has been computed [22] but we do not report it here. 
For both the universal functions corrections to the Gaussian value already at one-loop 
order have been found, at variance with Model A. 
In [22] the FDR has been computed, finding 

2 1 N + 8 



N ~ l ■ N ^- N K m 2-N)] 



(4-iV)(2-A0 ' 4(7V-1)2"-'- "/ijj +°( e2 )' ( 183 ) 

that, for TV = 1, is exactly the same as in Model A, X°° = 1/2(1 -e/ 12) + 0(e 2 ). This is 
really surprising, in particular looking at the complicate expression for general N. It is 
probably only a coincidence of one-loop computation, but only higher-loop calculations 
may clarify whether this is a deeper property or not. 



8. Model A dynamics of a weakly dilute Ising model 

A question of theoretical and experimental interest is whether and how the critical 
behaviour is altered by introducing in the systems a small amount of uncorrelated non- 
magnetic impurities leading to models with quenched disorder. 

The static critical behaviour of these systems with bond or site disorder (the effect of 
a random field is more complicate see, e.g., [157]) is well understood thank to the Harris 
criterion [158]. It states that the addition of impurities to a system which undergoes a 
second-order phase transition does not change the critical behaviour if the specific-heat 
critical exponent a p of the pure system is negative. If a p is positive, the transition is 
altered. 

For the very important class of the three-dimensional 0(A^)-vector models it is 
known that a p < for N >2 [67] and the critical behaviour is unchanged in the presence 
of weak quenched disorder (apart from large crossover effects, see, e.g., [159, 160]). 
Instead, the specific- heat exponent of the three-dimensional Ising model is positive [67]. 
Therefore the introduction of a small amount of non-magnetic impurities (dilution) leads 
to a new universality class (as confirmed by RG analyses, Monte Carlo simulations, and 
experiments, see [67, 161, 162] as comprehensive reviews on the subject) to which the 
weakly dilute Random Ising Model (RIM) belongs. The RIM is a lattice spin model 
with nearest-neighbour interaction Hamiltonian 

Tic = pipj SiSj , (184) 

(ij) 
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where Sj are the Ising spins at the site i, and pi are uncorrelated quenched random 
variables such that Pi = 1 with probability c (0 < c < 1 being the spin concentration in 
the lattice), p« = with probability 1 — c. Above the percolation threshold c*, the critical 
properties of the model (184) are predicted to be independent of the actual impurity 
concentration [67, 161, 162]. However to make this property apparent, the corrections 
to the scaling have to be properly taken into account in the analysis of numerical and 
experimental data (see, e.g., [163, 67, 159, 160]). 

The purely relaxational equilibrium dynamics (Model A) of this new universality 
class is qualitatively well understood [164, 165, 166, 167, 168, 169, 109]. The 
dynamic critical exponent z differs from the mean-field value already in the one-loop 
approximation [164], at variance with the pure model. This exponent is known up to 
two loops in a y/e [166] and up to three loops in fixed-dimension expansion [167] in 
three dimensions. Unfortunately the agreement between FT estimates (giving z ~ 2.2), 
numerical [169, 168, 61] (z ~ 2.6) and experimental [170] (z ~ 1.9 — 2.2) estimates is 
still quite poor. 

The non-equilibrium dynamics is instead less studied. Within the \/e expansion, 
6 was determined up to two-loop order [85] and the response function up to one 
loop, both for conserved and non-conserved order parameter [84]. A careful three- 
dimensional numerical simulation [61], taking properly into account the corrections 
to the scaling, showed that even is dilution independent and its value is (quite 
surprisingly, compared with z) in good agreement with the \fe expansion. In fact 
the obtained value X/z = 1.05(3) [61], using the scaling law X/z = d — 6', leads to 
& = 0.100(35) that has to be compared with the two-loop FT result ^ loops ~ 0.087 [85]. 
Also the autocorrelation scaling function obtained by means of numerical simulation 
[61] is in nice agreement with the RG prediction equation (92). The finite-size scaling 
in the non-equilibrium regime has been also investigated [171]. 

The time evolution of the weakly dilute Ising model is described by the stochastic 
Langevin equation (33) with 7i i— >• Ti^, where the static Landau- Ginzburg Hamiltonian 
Tiil>[<p\ has a space-dependent random temperature [172] 



^(Vy?) 2 + \{r G + ^(x))^ 2 + ^g oV 4 



(185) 



Here ^(x) is a spatially uncorrelated random field with Gaussian distribution 

s/Attw \_ 4w_ 

Dynamical correlation functions, generated by Langevin equation and averaged over the 
noise (, are given by equation (35) with Ti replaced by (we denote the resulting 
dynamical functional). 

In the analysis of static critical behaviour, the average over the quenched disorder ip 
is usually performed by means of the replica trick [173, 172]. If instead one is interested 
in dynamic processes it is simpler to perform directly the average at the beginning of 
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the calculation [174] 



P{ip) exp(-S^[p, <p\) = exp(-S[(p, <p\) 



obtaining the -^-independent action [85] (with v oc w) 

S[<p,0\= [d d x\ [ dt(p [dt<p + ^(r - A)tp - S10\ 



(187) 



+ 



dt (pip 3 - 



oo \ 2 

dt (pip 



;i88) 



3! J ~ rr 2 Vo 
Note that a non-local (in time) effective interaction term has been generated. 

The perturbative expansion is performed in terms of the two fourth-order couplings 
go and vq and using the propagators given in equations (112) and (113). The scaling of 
the response function is characterized by [21] 

A R = 1 \ — 

2 V 53 

and the correlation (see equation (87)) by [21] 



-7b + 0(e) , and F R (x) = 1 + 0(e), 



(189) 



£ =l + ^|(2-7-) +0 (e), 



CK J 2 V 53 



2 ( x) S T+x 



+ 0(e) . 



(190) 
(191) 



Note that, at variance with the pure model, the function Fq{x) gets a contribution 
already at one-loop order. Fr(v) remains 1 in apparent agreement with LSI although 
z 7^ 2, already at this order. 

The FDR for generic times is 



* q=0 (s/t) 



1 - 




6e 
53 



, 1 , I- st 
1 + - log -r 

2 s 1 + st 



+ 0(e) 



(192) 



X c 



+ 0(e), 



(193) 



Therefore, assuming even in this case the validity of equation (21), one finds 

1 1 l~6e 

2 ~ 4V 53 

that for e = 1 (that is the only physical relevant case, given that a p = in the two- 
dimensional Ising model) leads to X°° ~ 0.416. To this order it is not clear whether 
randomness really changes X°° in a sensible way or not. In any case this could not be 
safely stated from low-order computations since the y/e is known to be not well-behaved 
at d — 3 [67, 161]. However, a qualitative conclusion can be drawn from equation (192): 
X n= o has a quite strong dependence on the time ratio s/t, that, at variance with the 
pure model, should be easily identified in simulations. 

The relation between and in the RIM has been studied in more detail 

by Schehr and Paul [61]. Due to the quench disorder, the response function i? q has a 
power-law decay for q z t ^> 1 [84]. Therefore the argument leading to X^L Q = X™ Q does 
not apply straightforwardly in the form outlined in section 2.4. To clarify this point 
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Rq(t, s) and Cq(t,s) were calculated for generic momenta up to the first order in the 
y^-expansion. Then, via the Fourier transform, one can compute explicitly R x =o{t,s) 
and C x= o{t, s) and therefore X x= o(t, s). In doing that one faces the problem that i? q and 
C q [61] have a power-law decay for large q and fixed times, in contrast to the exponential 
one of the pure model. As a consequence one has to introduce an ultraviolet cut-off A , 
finding [61]: 

i? x=0 (t, S ) = K + ^log(t- S )](t- S r d / 2 Q , (194) 

C X=0 (M) = [A° c + ^log(* - 8 )]{t - sT +1 - d ' z [^j F c (s/t), (195) 

where A° RC are non-universal amplitudes (with A l R = A l c ). In particular A° RC explicitly 
depend on A and are not finite for A — > oo. Note that these expressions do not agree 
with the general RG predictions equations (91) and (92), because of the presence of 
the logarithmic term. However this fact is not surprising: Logarithmic corrections to 
the scaling behaviour are generically expected in the presence of quenched disorder, as 
shown by Cardy using conformal invariance in generic dimension [175]. 

We point out (as already noticed in [61]) that, up to one-loop order, the logarithms 
can be absorbed in the exponent of (t — s), giving rise to non-universal, cut-off dependent 
exponents a^^ = a + A l c R / A° c R . This possibility seems rather unlikely, but only a two- 
loop calculation can rule it out. We mention that the same kind of logarithmic behaviour 
has been found in pinned elastic interfaces near the depinning transition [176]. 

From equation (194) and (195) one easily derives the FDR in the real space [61] 

X x l = F x (t/s); with jr x ( u)= 2^±L + y|^y + 0( c ), (196) 

where Fx{u) interpolates between 1 in the quasi-equilibrium regime for u — > 1, and 
its asymptotic value for u — > oo given by lim^^ F x {u)~ l = X°° of equation (193). 
This result explicitly shows, at order y/e, that the asymptotic FDR for total and local 
magnetization are indeed equal, despite the power-law behaviour of i? q and C q for large 
q. This fact calls for a deeper understanding of the relation between X^L and X%L Q and 
of the fact that equation (21) seems to hold independently of the argument presented 
in section 2.4. 

9. Purely dissipative dynamics of a ip 3 Landau-Ginzburg theory 

So far we have only considered models whose static critical properties are described by 
the effective Landau-Ginzburg Hamiltonian given in equation (22). In all these 
cases 9 > and < X°° < 1/2. To understand whether these bounds on such 
universal quantities are general properties or they are accidentally connected with 
the diagrammatic properties of y? 4 theories, it is worth studying different effective 
Hamiltonians. In particular we consider a model whose static critical properties are 
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described by the <p 3 Landau-Ginzburg Hamiltonian 



HM = J d d x 



1 1 N 

' i,j,k=l 



(197) 



where the coefficients d^ specify the model and ip is a A^-component field. A variety 
of critical phenomena belongs to the universality class of equation (197). Examples 
are provided by the isotropic-to-nematic phase transition in liquid crystals [177], the 
A" + 1-state Potts model [178], its percolation transition (q — N + 1 — > 1) [179] and 
its limit q — > 0, representing electrical resistor networks [179]. The model with one 
component A^ = 1 and d m = i describes the Yang-Lee edge singularity [180]. Moreover 
a (/9 3 -like interaction is one of the most relevant in the effective field-theoretical action 
of the Edwards- Anderson spin glass [181]. 

The upper critical dimension for the universality class specified by equation (197) 
is d c — 6. Therefore it is possible to compute critical exponents and all other universal 
quantities in a e-expansion with e = 6 — d [182, 183, 184, 185]. Fixed-dimension results 
have also been provided [186]. 

However the field-theoretical study of the Hamiltonian (197) presents several 
potential pitfalls. For example in the e-expansion of the Potts-symmetric Hamiltonian 
a stable fixed point of the renormalized coupling constant g is found only for N < 7/3 
and its value g* diverges as A" approaches 7/3, whereas it takes imaginary values for 
N > 7/3. On the other hand the e expansion is not reliable for all values of A" for which 
g* is not sufficiently small [184]. This condition naively leads to the conclusion that the 
e-expansion makes sense only for A^ < 1. In spite of these problems the universality class 
of the percolation transition (N = 0) has to be considered free from pitfalls. This fact is 
corroborated by the remarkable agreement of field-theoretical estimates with numerical 
and experimental results [187] in d > 2 (in d = 2 it was claimed that the e expansion 
breaks down due to the relevance of tp 4 interactions [188]). 

The simplest dynamics that one can think of realizing for some of the model just 
described is the purely dissipative one, i.e., the Model A of section 3.4, defined by the 
Langevin equation (33). The dynamic critical properties can be worked out from the 
dynamical functional (35) with TL replaced by TC3. We point out that the dynamics just 
introduced is not the actual dynamics of isotropic percolation (reviewed, e.g., in [189]). 

The calculation of the one-loop non-equilibrium response and correlation function 
following a quench from the disordered state (introduced by means of Tio, as in 
section 3.5.2) to the critical point is straightforward and we do not provide any details 
but only the final results. Zero-momentum response and correlation functions satisfy 
the expected scaling laws (86) and (87) with exponents 

(a agrees with the known exponents [184, 190] through the scaling relation a = 
(2 — rj — z)j z) and scaling functions 

f «M = 1 - £ i(^4r +0 < e2 >- < 199) 
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F c (x) = 1 - e- 



'8(a-4/3) 
The resulting FDR is given by 

X°° = j[l+. 3a ~Le 



X-11 + 



X 



log(l-x 2 ) + C(e 2 ) 



+ 0(e 2 ). 



(200) 



(201) 



8(a-4/3) 

In the previous equations we introduced the three contractions 

dijkdjki = cxSii , diimdjmndkni = fidijk , dijkdku = 7<5y , (202) 

that completely characterize the model at one-loop level. Accordingly, for the models 
we are interested in, we do not report the coefficients d^ [178, 180], but only the 
contractions: 

Yang-Lee edge singularity: a = — 1, f3 — — 1, 7 = — 1, (203) 

N + 1 - state Potts Model: a = (N + 1) 2 (N - 1), (3 = (N + 1) 2 (N - 2), 7 = 0. (204) 

The obtained F R (x) disagrees at one-loop level with the prediction of Local Scale 
Invariance F R (x) = 1 [135]. This fact is not surprising, since at one-loop level, we 
have 2^2 (see the remarks in section 4.3). 

Let us discuss these results for some specific model. For the Yang-Lee edge 
singularity we found 







1- — e + 0(e 2 ) 
28 V ; 



— e + 0(e 2 ), X°° = - 
56 V ; ' 2 

i.e., > and X°° < 1/2. For the models with the symmetry of the q 
Potts model, instead 

1-N „ /0% 1 

2 



(205) 
N + 1-state 







-e + 0(e 2 ), X c 



3 1-JV , ' 

1 - -- — r^e + 0(e 2 ) 



7-3AT 



(206) 



8(7 -3A0 

For N < 1, where the perturbative expansion of such models is considered under control, 
we get > and X°° < 1/2. Conversely, for iV > 1 we apparently get X°° > 1/2 and 
< 0. However, such a result is not completely under control, given that one finds 
v < 1/2 [184] for TV > 1, in disagreement with the results for the q = N + 1-state Potts 
model [191]. 



10. Conclusions and Perspectives 

We reviewed the ageing properties of systems which are quenched from the high- 
temperature phase exactly to their critical points. It clearly emerges that the theoretical 
understanding of such phenomena is quite complete in comparison with the richer ageing 
effects in glassy phenomelogy. Nevertheless there are still some points that require 
further investigation in order to be completely understood: 

First, we mention the equivalence between the asymptotic FDR X°° defined in 
momentum and in real space, see equation (21). To our knowledge, the only general 
argument supporting this result is the one outlined in section 2.4 (originally proposed in 
[19]). However, a recent analysis of the dynamics of the weakly dilute Ising model [61], 
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suggest that probably such an equivalence has a wider range of applicability. It would 
be desirable to have a clear proof of it and to know the conditions under which one 
can (or can not) expect equation (21) to be valid, especially in connection the possible 
definition of the effective temperature (see section 2.4). 

Another interesting issue is the behaviour of these critical systems just at their 
lower critical dimensions (i.e., the dimension where the critical temperature vanishes, 
d — 1 for Ising-like systems and d — 2 for systems with continuous symmetries). 
Indeed, in passing, we mentioned [90, 94, 95] that different one-dimensional systems 
with Ising symmetry may have different non-equilibrium critical exponents (i.e., 6). 
In this context it has been recently argued [53] that a quench to zero-temperature at 
the lower critical dimension is not the dimensional continuation of a line of critical 
quenches in the (T, d) plane (as often implicitly assumed), but it is the continuation of a 
line of zero-temperature quenches, i.e., the system behaves as in the coarsening regime, 
although X°° 7^ (see for details [53]). Moreover the critical quench at the lower critical 
dimension has a very interesting aspect: Being the autocorrelation exponent a — djz + 1 
(see equation (91)) equal to zero, it is possible to write the FDR X x=0 as a function 
of the autocorrelation function C x= o as it happens in more complex instances of glassy 
behaviour (see, e.g., [7, 8]). 

Field-theoretical RG provides methods going well-beyond the expansion close to 
the upper critical dimension that has been used so far and that we have reviewed. It 
will be interesting to consider different approaches like expansions close to the lower 
critical dimensions, perturbative expansions in fixed-dimensions and non-perturbative 
approaches like the so-called Exact RG and 1/iV-expansion. These methods might 
lead to more accurate theoretical estimates for the universal quantities of interest and 
hopefully to a more general understanding of the phenomena. 

It is natural to wonder whether some of the results that have been reviewed 
carry over to glassy systems, which was the original motivation of the study of 
ageing behaviour in critical models. Recent numerical simulations of finite-dimensional 
id = 3, 4) spin-glass models exactly at the critical point [192] show that even in this case 
the response and correlation functions scale according to equations (91) and (92), as in 
critical non-disordered system, with a non-trivial value of X°° . On the other hand, for 
a quench in the low-temperature phase a completely different behaviour is expected [8]. 

Last but not least, we remark that the ultimate check of all the physical theories is 
provided by experiments. Unfortunately, to our knowledge, only quasi-one-dimensional 
systems with Ising symmetry have been considered so far [193]. Although very 
interesting, the physics of one-dimensional systems is too peculiar (d — 1 is the lower 
critical dimension of Ising systems) to provide a test of the theoretical results currently 
available. Experiments in two and three dimensions are required to check the relevant 
aspects of ageing dynamics at the critical point. In principle, this could be done by 
considering the same experimental systems as those considered when studying phase 
ordering — magnetic systems, nematic-liquid crystals, binary liquids, etc. [3] — exactly 
at their critical points. 
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